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Abstract We present a theory for the construction of renormalized kinetic 
equations to describe the dynamics of classical systems of particles in or out 
of equilibrium. A closed, self- consistent set of evolution equations is derived 
for the single-particle phase-space distribution function /, the correlation 
function C = {SfSf), the retarded and advanced density response functions 
T^fl,A _ Sf/S(f) to an external potential (f), and the associated memory func- 
tions S^'"^''^ . The basis of the theory is an effective action functional J7 
of external potentials (p that contains all information about the dynamical 
properties of the system. In particular, its functional derivatives generate 
successively the single-particle phase-space density / and all the correla- 
tion and density response functions, which arc coupled through an infinite 
hierarchy of evolution equations. Traditional renormalization techniques (in- 
volving Legendre transform and vertex functions) are then used to perform 
the closure of the hierarchy through memory functions. The latter satisfy 
functional equations that can be used to devise systematic approximations. 
The present formulation can be equally regarded as (i) a generalization to dy- 
namical problems of the density functional theory of fluids in equilibrium and 
(ii) as the classical mechanical counterpart of the theory of non-equilibrium 
Green's functions in quantum field theory. It unifies and encompasses previ- 
ous results for classical Hamiltonian systems with any initial conditions. For 
equilibrium states, the theory reduces to the equilibrium memory function 
approach used in the kinetic theory of fluids in thermal equilibrium. For non- 
equilibrium fluids, popular closures of the BBGKY hierarchy (e.g. Landau, 
Boltzmann, Lenard-Balescu-Guernsey) are simply recovered and we discuss 
the correspondence with the seminal approaches of Martin-Siggia-Rose and 
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o£ Rose. and we discuss the correspondence with the seminal approaches of 
Martin-Siggia-Rose and of Rose. 
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1 Introduction 

The present paper elaborates the foundations of a unifying approach for the 
construction of renormahzed kinetic equations of classical systems of par- 
ticles in and out of equilibrium. The statistical dynamics is formulated in 
terms of self-consistently determined single-particle phase-space distribution 
function, density response, correlation and memory functions. The memory 
functions satisfy a functional equation that can be used to devise systematic 
renormahzed perturbative approximations. Important topical applications in- 
clude the systematic generalization of the Boltzmann equation to dense gases 
and liquids, the kinetic theory of strongly coupled plasmas and Hamiltonian 
systems in general, the proper macroscopic description of fluids beyond the 
Navier-Stokes equations, and the physics of plasma turbulence. 

Before embarking on the detailed description of the theory, we recall a 
few of the most striking developments and unresolved issues in kinetic theory 
that motivated this work. 



1.1 Background. 

Kinetic theory attempts to carry out a complete statistical description of the 
macroscopic dynamics of fluids in terms of the underlying microscopic inter- 
actions of its constituents. The foundation of modern kinetic theory dates 
back to 1872 when Boltzmann published his famous equation for the single- 
particle distribution function for dilute gases [HHIII]. The solution of the 
Boltzmann equation, derived independently by Chapman and Enskog around 
1915, achieved two basic goals of kinetic theory. That is, it (i) established the 
connection between the microscopic dynamics and the macroscopic equations 
of hydrodynamics and (ii) provided explicit expressions for the transport co- 
efficients in terms of the molecular parameters. Unfortunately, since it is 
based on the assumption of uncorrelated binary collisions, the Boltzmann 
equation gives a satisfactory description of the behavior of sufficiently low 
density gases with short-range interactions only, and is inadequate to describe 
dense gases, liquids and plasmas. 

The generalization of the Boltzmann equation to more dense fluids and 
to plasmas confronts the full complexity of a highly correlated many-body 
problem. This became very clear when, in the early 1960's [4], it was real- 
ized that all the remarkable systematic developments made until then, which 
were based on some form of perturbation theory of the Liouville equation 
or the BBGKY hierarchy, were plagued by unphysical divergences [5j|6l[7]. 
Those findings, together with others permitted by the advent of computer 
simulations, shifted attention towards the development of a renormalization 
(regularization) procedure, which would encompass the "theoretical" diver- 
gences and take proper account of the correlated collisions, whereby particles 
interact via the effective (renormahzed) potential that integrates the aver- 
age effects of the medium on the bare interactions [SUH] ■ Additionaly, such a 
theory should yet retain the desirable properties of the Boltzmann equation 
and be tractable for practical calculations. 
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Such a renormalized kinetic theory was first presented by Mazenko in 1974 
using the memory function approach [12] . His so-caUed "fully renormalized" 
theory apphes only to fluids in or near thermal equilibrium by means of the 
fluctuation-dissipation theorem. At equilibrium, the distribution function / 
is known (the Maxwcll-Boltzmann distribution) and the statistical dynamics 
is encoded in the equilibrium two-point time-correlation function Ceq(l,l') 
of the microscopic phase-space density. It was shown in the 1960's [IS] that 
an exact evolution equation for Ceq can be written in the form of a general- 
ized Langevin equation in terms of a memory function kernel S^q- Mazenko 
devised an iterative scheme for its calculation that can be used to develop 
systematic renormalized approximations of Seq- 

The development of a renormalization procedure that can describe non- 
equilibrium states faces additional challenges [13 . A first major development 
in that direction was done in 1973 by Martin, Siggia and Rose (MSR) [T4] . 
MSR elaborated a procedure to derive self-consistent approximations for cal- 
culating the statistical dynamical properties of a classical random variable 
or field whose time dependence is governed by a nonlinear differential evo- 
lution equation. In other words, they devised a systematic approach to the 
ubiquitous statistical closure problem, which occurs as a consequence of the 
statistical averages and the non-linearity of the fundamental evolution equa- 
tion. Noting that satisfactory approaches to the closure problem had been 
successfully developed in quantum field theory, MSR recast the classical prob- 
lem into a quantum-mechanical-looking problem and borrowed techniques of 
quantum field theory to derive a renormalized perturbation theory. To ac- 
complish the recast, MSR introduced an additional field operator that does 
not commute with the fundamental field and allows the calculation of both 
correlation functions and response functions simultaneously, where the latter 
describe the response of the fluid to the injection of particles. 

Although the basic problem in kinetic theory also amounts to a closure 
problem, it was then recognized that its specificities did not fit into the 
framework of MSR theory; in particular, its natural field variable, the micro- 
scopic single-particle phase-space density, is discrete, singular and its statis- 
tics is strongly non-Gaussian. In 1979, Rose [TS] developed another approach 
specifically for dealing with kinetic theory. Like MSR, Rose also recast the 
basic problem in a quantum-mechanical form but this time using an occu- 
pation number representation in classical phase-space |16j . The response to 
the injection of particles included in MSR is naturally incorporated here. 
The implications of Rose's paper to kinetic theory have not yet been fully 
explored. In addition. Rose points out that the steps necessary to reduce his 
general formalism to Mazenko's theory of fluids in thermal equilibrium are 
not apparent. 



Finally, very recently, Mazenko [ITj published a new fundamental "the- 
ory of statistical particle dynamics that solves the chronic problem of self- 
consistency" . In that paper, the theory is presented for the Smoluchowski 
dynamics and the author announces that its extension to Newtonian dynam- 
ics is under way. From the information available in |17) . it is likely that the 
approach presented here is still different from Mazenko's new theory. 
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1.2 Present work. 

In this paper, we present an alternative, arguably simpler and elegant ap- 
proach to the derivation of renormalized kinetic equations for classical fluids 
that applies to both equilibrium and non-equilibrium states and encompasses 
the previous works on this topic. For simplicity of the exposition, we shall 
consider a system of identical point particles mutually interacting via a pair- 
wise additive central potential; however, no peculiar assumption is made on 
its range and strength, so the theory can be applied to investigate systems 
ranging from neutral liquids with short-range interactions to classical plasmas 
with the Coulomb interaction. Moreover, the theory is derived independently 
of the nature of initial state of the system, and the influence of initial cor- 
relations on the system dynamics, usually neglected or difficult to include in 
the previous works, is treated explicitly. 

A closed, self-consistent set of evolution equations for the phase-space 
density f{r,p,t), the two-point correlation function C, and the retarded and 
advanced response functions x^'^ is derived in terms of time history integrals 
that involve three memory functions X'^'^''^. Loosely speaking, the correla- 
tion and response functions describe the dynamics of emission and absorp- 
tion of phase-space density fluctuations in the fluid, which in turn determine 
the effect of collisions on the dynamics of the distribution function /. The 
memory response functions describe how many-body effects affect, i.e. renor- 
malize, the time-dependent propagation and lifetime of density flucutations, 
and as such play a role similar to self-energies in quantum field-theory. 

To obtain this closure, we construct an action functional of external 
potentials cj) that contains all information about the dynamical properties of 
the fluid. In particular, its derivatives with respect to (f> successively generate 
the phase-space density / and all the correlation and response functions. The 
latter are coupled through an infinite hierarchy of evolution equations, which 
subsumes the traditional BBGKY equations. The hierarchy is formally closed 
using traditional closure (renormalization) techniques involving the Legendre 
transform r[f] = —f2[(j)] + J (j>f of the effective action 

The basis of the present approach is the extension to classical systems 
of the closed-time contour idea originally introduced in 1961 by Schwinger 
|18j in his seminal paper on the Brownian motion of a quantum particle. 
Schwinger's idea was later fruitfully used to develop the theory of non- 
equilibrium Green's functions in quantum field theory |191l20j . Two fully 
equivalent formulations of the theory are presented. In the first formulation, 
the theory is expressed in terms of time-dependent quantities taking their 
values on the closed-time contour. The effective action J7[0] is a functional of 
external potentials (j> that linearly couples along the closed-time contour to 
the single-particle phase-space density. The first functional derivative with 
respect to 4> generates the phase-space distribution function / = Sf2/S(j) and 
the second derivative S'^n/Scj)^ equals a quantity that combines both the two- 
point correlation function C and the density response functions x^'"^- In the 
second formulation, the theory is recast into a form that directly generates 
the physical quantities of interest. The effective action is written as a func- 
tional of two external potentials, fl[4>p, (I)a\, of the physical time variable t 



6 



and has the remarkable derivatives / = , C = s/^2t>^ ' ~ s^Jh ^"^^ 
A _ s'^n 

The main purpose of this paper is to lay down the general foundations 
of the theory without reference to a specific physical system or model. Its 
application to physical problems would require substantial additional work, 
some of it being under way. The paper is organized as follows. In Sec. [2] the 
basic definitions and notations are introduced. To provide a self-contained 
presentation of the work, we recall several basic results of classical statistical 
mechanics that are used thereafter. The original contributions of this work 
really begin into Sec. [3] where we develop the closed-time formulation of 
our theory and continue in Sec. |4] and [5] where we recast the results into 
equations that involve quantities depending on the physical time. For clarity, 
the details of the proofs of the calculations are given in the appendices. 
Throughout these sections, we discuss the implications and merits of the new 
approach presented. We study the special case of a fluid at equilibrium, make 
contact with Mazenko's renormalized theory, show the correspondence with 
popular kinetic equations (Landau, Boltzmann, Lenard-Balescu-Guernsey), 
and study the correspondence with MSR and Rose's theories. Also the present 
approach naturally shed light on analogies and differences between quantum 
and classical kinetic theory, which we highlight throughout the paper and in 
appendices. The comparison allows us (i) to comprehend the longer resistance 
of classical statistical mechanics to renormalization, and (ii) to justify some of 
the ingenious approaches used by previous works like MSR. For convenience, 
a summary of the main components of the theory is given in SecjGj 



2 Basic definitions 

In this section, we introduce the physical systems under consideration. For 
completeness, we also recall several basic notions of statistical dynamics used 
thereafter such as the correlation functions, the linear response functions to 
external perturbations, and the closure problem. Most of the material is 
discussed in textbooks, e.g. (SHU], and is recalled here for completeness. The 
reader already familiar with these notions may wish to skip ahead to SecjS] 
where the original contribution of this paper begins. 



2.1 Physical system 

We consider a system consisting of N identical point-particles of mass m 
enclosed in a volume V of the d-dimensional space M''. We assume that the 
system dynamics is governed by the laws of classical mechanics under the 
total, possibly time-dependent Hamiltonian Htot- Thus, the position rj{t) 
and momentum pj (i) of the j-th particle at a time t evolve according to the 
Hamilton equations, 

drj _ dHtot(x,t) dp J _ dHtot{x,t) 
dt dpj ' dt dvj 
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where a; = (ri, . . . , r^r; pi, . . . , pjv) e x M''^. These equations ([T]) are to 
be solved subject of 2dN initial conditions on the coordinates and momenta 
x{to) = xq at an initial time to- 

We assume that the initial state xq of the system is imperfectly known but 
can be characterized by a distribution Fo{xq) defined such that Fo{xo)dx() is 
the probability that the system is initially in a microscopic state represented 
by a phase-space point of volume dxo around xq; in the following, _fo is as- 
sumed to be normalized, i.e. f dxQFo{xo) = 1. Most importantly, the present 
theory is derived independently of the nature of initial state Fq of the system. 

Thus we are interested in the dynamics at time t > to of a statistical 
ensemble of independent systems, each of which is a replica of the system 
defined above, and initially distributed according to Fq; the dynamics prior 
to the initial time to is not of interest to us. The uncertainty on the initial 
conditions will be the only element of statistics in the present theory. Once the 
non-equilibrium initial state is specified, the time-evolution is deterministic 
and completely determined by the Hamiltonian. 

As time increases, the initial distribution of states Fq evolves into the 
distribution F{x,t) at i > tg according to Liouville equation, 

^^F{x,t) = [Htot,F]pg{x,t), (2) 
where [A,B]pb denotes the A^-particle Poisson bracket 
[A B]pB = Y (— ,dB_ _d^ dB\ 

We assume that the total Hamiltonian can be written as 

Htot{x,t)^HN{x) + H,,t{x,t). (3) 

The time-independent term iJ^r characterizes the system in the absence of 
external, time-dependent perturbations, and is of the form, 

N I ^ 

= ^/io(rj,p,) + - E v{\r,-r,\). (4) 

j = l i^j=i 

The particles interact according to a pair-wise additive central potential 
v[r) and the whole system is possibly confined by a static potential wo(r)- 
Throughout the paper, no peculiar assumption is made on the range and 
strength of the interaction potential v other than being differentiable for 
r > 0. In particular, the theory applies to both short-range potentials as 
encountered in neutral liquids [12] and also to the long-range, Coulomb po- 
tential of importance to plasma physics; in the latter case, a background 
potential may be included in vo to ensure electrical neutrality [23]. 
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The term H^^t in Eq.Q represents the possible interaction to an external, 
time-dependent perturbation (f>o(r,p,t), 

Hext{x:t) ^^<j>o{rj,Pj,t) ^ // drdpN{r,p,t)(t)o(r,p,t) , (5) 
3 = 1 '''' 

which linearly couples to the phase-space density iV(r,p,i:) (see Eq.([7|) be- 
low.) 

The general definitions given above encompass a large class of systems, in- 
cluding for instance (i) the dynamics of a fluid at equilibrium when H^xt = 
and Fq — feq is an equilibrium (Gibbs) distribution function, (ii) the relax- 
ation dynamics to equilibrium when Fq ^ /eg and H^xt = 0, (iii) the out- 
of-equilibrium dynamics of a fluid in the presence of time-dependent fields 
when Hext 7^ 0. 



2.2 Closure problem, correlation and response functions 

2.2.1 Fundamental field variable and correlation functions 

A dynamical variable is a quantity A that depends parametrically on the 
particles' trajectory x(t), and therefore on the initial condition xq. As a 
consequence of the statistical description, A no longer has a definite value 
but is instead characterized by its average over all possible initial conditions 
weighted by Fq, 

{A){y)^ J dxoFo{xo)A{xo;y); (6) 

here y denotes the possible dependence on other parameters. 

A fundamental dynamical variable in classical kinetic theory is the mi- 
croscopic single-particle phase-space density 

N 

7V(r, p,t) = J2Hr- r, (t)) S (p - p, (i)) , (7) 
j=i 

where (r, p) ^ V x W^. The density N{r,p,t) plays a role similar to a field 
in field theory or of a random variable in the theory of stochastic processes; 
here N{r, p, t) is random in that its value varies due to its sensitivity to ini- 
tial conditions. The ensemble average of A^(r, p, t) defines the single-particle 
distribution function 

/(r,p,t) = (iV(r,p,t)), (8) 

and has the meaning of the probability distribution of particles at position r 
with momentum p at time t. For an unconfined fluid {vq — 0) in thermal equi- 
librium at temperature T, f is independent of space and time, and reduces 
to the Maxwell-Boltzmann distribution /(p) = ne'^'/^"'''^'^ /(2TrmkBTf^'^ 
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where n — N/V is the particle density. In non-equihbrium fluids, / generally 
depends on both r and t. 

As is typical in a statistical theory of a random variable N{1), an im- 
portant role in the theory is played by the fluctuations around its averaged 
value, 

dN{l)=N(l)-{N{l))^N{l)-f{l), (9) 

and by the associated n-point correlation functions (or cumulants) 

C(") (1, . . . , n) = {SN{1) . . . SN{n)) . (10) 

For the two-point correlation function, we write C(l,2) = C'-^^(l,2). Equal- 
time correlation functions, i.e. ti — ■ ■ ■ — tn, will be denoted by a bar as 
C(")(a;i,...,x„;t). 

As was emphasized by Rose [TS], a statistical theory based on a dis- 
crete (i.e., a sum of delta functions) random variable such as N{1) presents 
a serious technical difficulty in that it is strongly non-Gaussian, even for 
non-interacting particles, since the equal-time cumulants C'"' are all non- 
vanishing and singular. Indeed, in terms of the conventional n-particle cor- 
relation Qn introduced in kinetic theory from the cluster expansion of the 
reduced distribution functions (see appendix |A|, we have 

C^^\X,X' ,t) ^ 5{X - X')f{X,t) + g^{X,X';t) (11) 

and 

(7(3) (X, X', X"- 1) = 5{X - X')5{X' - X")f{X, t) 

+ [5{X - X')g2{X', X"; t) + cyclic permutations] 
+ g3iX,X',X";t) (12) 

and so on, with X = (r,p). Thus, even in the absence of three-particle 
correlations (173 ~ 0), the three-point cumulant is non- vanishing because of 
particle self-correlations, and this remains true at all times t > tg. This 
result precludes the applicability the Wick theorem so useful in conventional 
statistical field theories. 

2.2.2 Field equation and closure problem 

For convenience we introduce the shorthand notation in which a field point 
and the time variable are designated by a single number, i.e. n = (r„, p„, tn) 
(not to be confused with the particle labels n); for n = 1, we often drop the 
subscript so that 1 = {r,p,t). The delta function (5(1 — 2) denotes ^'^^^(ri — 

r2)5('HPi-P2)5(ii-i2). 

From the Hamilton equations ([T]), it is straigthforward to show that the 
phase-space density evolves according to 



— N{l)^LiN{l)+ J d2Li2iV(l)iV(2). 
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where the summation is defined as 



dl 



dri / dpi 



dti 



Here, 



dhjl) d* 
dvi 9pi 



dh{l) d* 
9pi dvi 



{h,.){l) 



(13) 



is the single particle part of the Liouville operator with the total single- 
particle Hamiltonian = /iQ + ^io- In the second line of Eq.( 13 1, it is expressed 
in terms of the single-particle Poisson bracket 



{a,6}(r,p) = 



da db da db 
dr dp dp dr 



(14) 



for any function a(r, p) and 6(r, p) in R^'^ (note the curly brackets notation to 
avoid confusion with the A^-particle Poisson bracket [•, ■]pb defined earlier.) 
Finally, 



il2» 



d , 



d» 

dpi 



d* 

dp2 



(15) 



is the interaction part of the Liouville operator with Vi2 ~ u(ri — r2)(5(ii — ^2)- 
Those notations are in widespread use in the literature on classical kinetic 
theory. In the theory presented below, however, it will be more convenient 
to rewrite Eq.(13) as 



d_ 
dfi 



Li Nil) 



d373(l,2,3)7V(2)iV(3) 



^73(l,2,3)iV(2)iV(3) 



in terms of the "bare interaction vertex' 

d 



dp I 

d . 

OTi 



(wi2^(l-3)+wi3<5(l-2)) 



(16) 



(17) 



73(1,2,3) 
W12 

which is symmetric in its last two arguments. The second line of Eq. ( 16 1 illus- 
trates the summation convention over repeated indices of dummy variables 
that we shall use throughout the paper. 

Averaging Eq.(16l over the initial conditions, we obtain the equation of 
evolution of the phase-space distribution /(I), 



jd_ 
dh 



Li 



/(l)-{«"/(l),/(l)}=-73(l,2,3)C(2,3). 



(18) 
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where u"''^ is the mean-field (a.k.a. Vlasov) potential 



(r,i) = / d2/(2){)(l - 2) 



(19) 



Equation (18 1 couples /(I) to the (equal-time) two-point correlation func- 
tion C(l,2j7and corresponds to the first equation of the ordinary BBGKY 
hierarchy; the right-hand side (rhs) is usually referred to as the collision in- 
tegral. Si mila rly, as a consequence of the quadratic non-linearity of the field 
equation (16 1, the evolution of C^"^ at each order n involves the next higher- 
order correlation (7*^"+^^; those equations can be straightforwardly derived 
using the equation for the density fiuctuations 6N obtained by subtracting 
Eqs.(16) and (18). Thus, for the two-point correlation function, we find 



_d_ 



Li 



where 



C(l, 1') - ^™/(l, 2)C(2, 1') = ^73(1, 2, 3)C(3)(2, 3, l') 



i:"/(l,l')=73(l,2,3)/(2)<5(l'-3) 



(20) 



The hierarchy of evolution equations for the C'"^ is different from the or- 
dinary BBGKY hierarchy, which instead involves the correlation functions gn 
(a quick summary of BBGKY is given in appendix |A)) The BBGKY can be 
recovered by considering the equal-time limit of the equations for the C*-"^'s. 
The "simplest" closures of the BBGKY hierarchy lead to the most famous ki- 
netic equations, in particular the Vlasov equation, the Landau equation [10] 
for so-called weakly-interacting gases, the Boltzmann equation [T] for dilutes 
gases, and the Lenard-Balescu-Guernsey equation [11] for weakly- coupled 
plasmas. Those closures rely, among additional hypothesis, on a systematic 
ordering of the correlations gn in terms of an adequately chosen small pa- 
rameter [B]. As summarized in tab le [ij each closure amounts to keeping only 
certain contributions in Eqs.( 11 )-( 12 1. In particular, these popular closures 
systematically neglect the n-particle correlation functions gn with n > 3, 
which is inappropriate whenever correlations are strong. A major challenge 
of kinetic theory is to incorporate the effects of those neglected terms in a self- 
consistent manner. In Sec l4.5[ we shall see how the present theory addresses 
the problem. 



2.2.3 Linear response functions and fluctuation-dissipation theorem 

We now consider the linear response of the system defined above to weak, 
external perturbations (j) that couples linearly to the single-particle density 
N{X,t) as 



Hpert{r,p,t) ^ / drdpN{r,p,t)6(j){r,p,t) . 
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Closure 


(5(2) 


C(3) 


Vlasov 








Landau 


5{x - x')f{x,t) 





Boltzmann 


5{x - x')f{x,t) 


[S{x - x')g2{x' ,x"-t) + c.p.] 
+S{x-x')5{x' -x")f{x,t) 


LBG 


g2{x,x';t) + 5{x - x')f{x,t) 





Book-Frieman |47j 


g2{x, x'; t) + S{x — x')f{x, t) 


[S{x~x')g2{x',x"-t) + c.p.] 
+S{x~y)Siy~z)fix,t) 



Table 1 The most popular closures o f cla s sica l statistical physics amount, among 
additional hypothesis, to replace Eqs. |ll| -( |l2| for C''^''''^' by the terms displayed 
in the table. Those approximations systematically neglect the n-particle correlation 
functions g„ with n > 3, which is inappropriate whenever correlations are strong. 
One of the challenge of kinetic theory is to incorporate the effects of those neglected 
terms in a self-consistent manner. 

For comparison, the term g2 kept in the Lenard-Balescu-Guernsey (LBG) equation 
but discarded in the Boltzman equation allows the renormalization of the bare in- 
teraction into a dynamically screened two-particle interaction between the particles 
of the plasma. On the contrary, the contribution to C^^' kept in the Boltzmann 
approximation but discarded in the Lenard-Balescu-Guernsey equation, is respon- 
sible for the bare two-particle interactions (large-angle scattering) describing the 
two-body collision in a dilute gas. 



By linearizing with respect to S<j) the Liouville equation Q with Hfot replaced 
by Htot + Hpert , we find that to first order in the external perturbation the 
phase-space distribution /(r. p, t) is modified from its value in the absence 
of perturbation by the amount, 



5f{l) = J dl'x^(l;l')W), 



in term of the retarded response function, 

X'^ih 1') = ^ - eit - mt' to){[Nil),Nil')] ^^), (21) 

for all t,t' > tQ. The retarded response function is not symmetrical in its ar- 
guments but instead x^(l'i 1) = X"^(l: where is the advanced response 
functions, 

X^(l; 1') = -e{t' - mt - to){[Nil),N{l')] . (22) 

In the Poisson bracket, N{1) is evaluated along the unperturbed trajectory 
governed by Htot (interaction representation). Other useful properties in the 
derivation of the present theory are given in appendix [Bj 

The two-point correlation C and response functions x^'^ play a funda- 
mental role in the present theory. The correlation function gives information 
about the likehhood of fluctuations in the fluid and determines the effect 



of particle collisions in the collision integral of Eq.(18l The response func- 



tions carry the dynamical information on how fluctuations are propagated 



13 



in time. In classical mechanics, those functions are related to averages of 
dynamical variables of different nature, namely the averaged product of field 
variables for C and their averaged Poisson bracket for x^'^. This result is to 
be contrasted with the equivalent result in quantum mechanics where both 
correlation and response functions are related to the average of products of 
the same dynamical variables, on which more is said later. 

Only for the special case of a system in thermal equilibrium [Fo = feq and 
Hext = 0), the response and correlation function are simply related according 
to 

x''{i,i')^mt~mt'-to)-^c{i,r), 

with /? — l/ksT. This result is known as the fluctuation-dissipation theorem 



3 Renormalization on the closed-time contour 

The present theory is based on the definition of a generating functional 
that contains all the information about the dynamics under investigation. 
Here 4> represents some external time-dependent potential, which perturbs 
the system and, just like (j)o in Eq.([5]), couples linearly to the microscopic 
phase-space density in the total Hamiltonian. The functional is designed such 
that its derivatives with respect to (f) generate the phase-space distribution 
/, the correlation and response functions C and x''^^^ of the system. 

The functional f^[4>] plays a role similar to the grand-potential 12 [wq] in 
the density functional theory of classical fluids used to study their static 
properties |22j : in this case, the first derivative with respect to the confin- 
ing potential vq gives the particle density 5n[vo]/6vo{r) — p(r) and the sec- 
ond derivative equals the density correlation function Sf2[vo]/Svo{r)Svo{r') = 
{5p{r)Sp{r'))eq [22]. The existence of such a generating functional for the 
calculation of dynamical properties is far from obvious. For instance, if one 
imagines a straightforward extension of the static theory in which i7[0o(l)] 
is a functional of time-dependent external potentials 4)o{r,p,t) that gener- 
ates the phase-space density Sil/S(j)o{l) — /(I), one is immediately led to a 
contradiction. Indeed, the symmetry of the second-order derivative requires 

5 5f2 5 5Q 



<5</>o(2) 50o(l) <50o(l) ^0o(2) ' 
while causality requires 

^ ^ _ Mil _ 2. , n,^ 

^</.o(2)<5</,o(l) -J<^o(2) ^'''^^^ ^'''^ 

since, say when <i > t2, X^(l,2) = {[N [1) , N (2)] p g) and x^(2, 1) = (see 



Eq.(21|). Both requirements could never be satisfied together and therefore 



the search for such a functional is doomed to fail. 
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In order to cope with the problem, we shall define f2 over a wider set 
of dynamics using the idea of closed-time contour originally introduced by 
Schwinger in the early 1960's to deal with quantum systems [T^. Because 
quantum mechanics propagates probability amplitudes instead of probability 
densities, the mere notion of closed-time contour arises more manifestly in 
quantum than in classical mechanics. Thus, if at time t > to a quantum 
system is in the state \'I'{t)), the expectation value of an observable O at 
that time can be expressed as 

(o)(i) = {m\o\m) (23) 

- {^{to)\U{to,t)dU{t,to)mto)) (24) 
in terms of the propagator U{t,t') from time t' to t, and the initial state 



|tf'(io)) at time to - If we read the time arguments of the propagators in Eq.( 24) 
from right to left we may say that the system evolves from Iq to time t after 
which the operator O acts; then the system evolves back along the real axis 
from time t to Iq. Schwinger first noticed that one could imagine that the 
forward and backward time evolutions are governed by different dynamics 
(e.g., if a different external potential acts on each branch), and showed how 
this can be turned into an ingenious device to probe (generate) the properties 
of the physical system under investigation simply by suitably perturbing its 
dynamics along the closed-time contour that goes from tg to t back to tg. 
The closed-time contour approach was further developed by Keldysh ^25j and 
many others to obtain an elegant and powerful treatment of non-equilibrium 
quantum systems in terms of non-equilibrium Green's functions p!^[^E5] . 



In contrast, in classical mechanics, the expectation value (23) becomes 
{0){t) = f dxoF{xo)0{x{t)) , 



and the concept of closed-time contour is much less apparent. Nevertheless, 
we shall show how a similar idea can be fruitfully developed to the classical 
case. 

First we shall define the notions of closed-time contour and extended 
dynamics. Then we shall use these concepts to define the action functional 
and investigate its generating properties. Finally we shall show how 
J7[(/)] provides a formal solution to the closure problem. 



3.1 Closed time contour 

We shall extend the dynamics under investigation along a closed path in 
time that, as illustrated in Fig.([T]), goes from the initial time to to some time 
<Af in the future and back to the initial time tg. The value of Im can be 
chosen arbitrarily as long as it is larger than the largest time of interest; for 
definiteness, we shall set tM = +oo in the following. The closed-time contour 
can be conveniently defined by parameterizing the physical time t = t^r) in 
terms of a time variable r in such way that if r runs from an initial to 
a final r/ then t(r) monotonically increases from the initial time to to the 
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^time 



Fig. 1 Schematic illustration of the closed-time contour. For any time t G [to,tM], 
r+ is the time on the forward branch c+ and t_ is the time on the backward branch 
c_ of the closed-time contour such that t = t(r+) = t{T-). 



i(T) 




Fig. 2 Illustration of the extended dynamics x{t) on the closed-time contour. For 
a physical potential, the dynamics retraces the forward trajectory on the return 
branch. 

maximum time at tm and then monotonically decreases from Im back to 
io [22] ■ In fact, the chosen parametrization is inconsequential since the final 
results are independent of it. 

In the following, for any physical time t G [to,tM], we will denote by 
the unique time on the forward branch and r_ the unique time on the 
backward branch of the closed-time contour such that t = t(T+) = i(T_). 
We will also distinguish certain quantities defined on the closed-time contour 
variable r by a tilde. Thus we write 1 — (ri,pi,Ti) — {Xi,ti) or simply 
1 = {X,t) when no confusion is possible. In these notations, d{l — 2) = 
5(3)(ri-r2),5(3)(pi-p2)5(ri-r2). 

3.2 Dynamics on the closed time contour 

Let us first imagine that the dynamics x{t) — {rj(t), Pj(t)}_^._^ ^ under in- 
vestigation over the physical time interval to < t < i^/, or equivalently 
Ti < T < Tm, extends to the backward branch tm < t < of the contour in 
such way that, as illustrated in Fig.(l2]), the extended dynamics retraces back- 
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ward the same phase-space trajectory along [tm, t/]. In other words, if tjir) 
and Pjir) denote the positions and momenta of this extended dynamics, we 
require that for all r € [Ti,ry], 

f,(r)=r,(t(r)) , p,(r) = p,(t(r)) . (25) 



By differentiation of Eq.(25l with respect to t, we easily find that the ex- 



tended dynamics is Newtonian and governed by the Hamiltonian, 



N 



n^„{x,T) = t'iT)Htoti£,tiT)) ^ t'ir) ^HN{S:) + Y,^Mrj.Pj,tiT))j (26) 

with X — (f 1, . . . , f^r; pi, . . . , ptv). This extended dynamics then can be re- 
garded as a special case of a larger set of dynamics on the closed-time contour 
and characterized by the Hamiltonians 

H^{x,T)=t\T) |^i?jv(i) + ^ </-(?„ P„r)^ . (27) 

= i'(r) (^Hn{x) + j) dxdTU{x, T)(j){x, r)^ 

where 4> is any external, time-dependent potential defined on the closed-time 
contour, and 



N 



A/-(r, p, r) = 5 (r - f , (r)) 5 (p - p, (r)) 



is the single-particle phase-space density. Just like in calculus it is often 
very useful to consider a real function of the real variable as a complex 
function of the complex variable, considering the physical dynamics under 
consideration as an element of the set of extended dynamics will allow us 
to derive properties of the physical dynamics that are difficult to obtain 
otherwise. 

We can distinguish between two types of external potentials </>, namely 
physical and non-physical potentials; we shall denote physical potentials with 
a subscript p, e.g. (pp. A physical potential is identical on both branches of 
the contour, i.e. (/)p(r+) = 4'p{T-), while a non-physical potential is not. Any 
potential with real, physical origin such as (jiQ in Eq.([5]), defines a physical 
potential through 4>p{t) — (j)o{t{T)). By construction, under the influence of 
a physical potential, the extended dynamics retraces its forward trajectory 
backwards on the return branch of the closed-time contour, while for a non- 
physical potential, forward and backward trajectories are generally different. 
Both cases are illustrated in Fig. [2] 
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3.3 Generating action functional 

The extended dynamics governed by can also be described by the La- 
grangian 

C^{f,v,T) =p-v-n4,{f,p,T), (28) 

with the velocity v — df/dr and the short-hand notation f = (fi, . . . , rAr). 
Given an initial condition xq and the corresponding trajectory (f(r), 'D(r)), 
we define the total action over the closed-time contour as 

/•^/ 
dr C^{f{T),v{T),T) ~Sb[(j),xo], 

-\-oo 

dTn{T)C^{f{T),i{T),T) -5f,[0,Xo], (29) 

-oo 

with the boundary term, 

Sb[(j),Xo\ = ^ {p[Ti) + p[Tf)) {r{Ti) - r{Tj)) . 

The role of Sb is to cancel boundary terms that arise when performing the 
variations of the first term with respect to (j). In the second line, we include 
the integration range [r^, Tf] in the integrand through the "window" function 

7T(t) = e{T - n)B{Tf ~t)=9 (i(r) - to) (30) 

where 6 is the (Heaviside) step function so that i7(T) = 1 when < 
T < Tf and n{T) = when r < or t > Tf. As we shall see, this appar- 
ently insignificant rewriting in fact allows us to self-consitently include the 
contribution of initial conditions within the theory. 

We then define the effective action functional Sl[(j>] such as 

e-i'^I'^l = /eJ^I^'-ol) =Z[(/)]. (31) 



Here s is a fixed parameter with the dimension of an action, or energy times 
time, to make the exponent dimensionless. Its value may be set to the quan- 
tum of action h if one regards e^'^I'^'^o] a,s the classical contribution to the 
total quantum amplitude to go from points i(ri) = xq to x{Tf) [28j . However 



the identification s — h is not required since Eq.(31) could also be defined 



without reference to its quantum origin. Accordingly, we find that the actual 



value of s is inconsequential since all the general results derived from Eq.(31 ) 
are independent of it. (The situation here is reminiscent of the problem en- 
countered in the early days of classical statistical mechanics [55], where an 
arbitrarily chosen size of the unit cell in phase-space was introduced in 
order to count the accessible classical states; while the physical laws were 
not affected by the value of /iq, the arbitrariness was later removed using the 
principles of quantum mechanics, leading to /iq = fi.) It may for instance be 
useful to chose s to be purely imaginary, s — is , and write 



= -slnZ[0] = -sln(e 
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in order to circumvent difficulties related to complex logarithms, to ensure 
nice properties of the effective action (e.g. convexity), etc.; we shall not delve 
into these technical points here. 

We remark that the action functional vanishes at any physical potential 

fi[<i)p] = , V0p . 

Indeed, the total action S vanishes 

S[(j)p\X(3\ = , Vxo , 

since the contributions over both branches of the closed-time contour have 
opposite sign by construction, 

j-TM j-Tj 

/ dr C^{t) ^ - I dr C^{t) , 

J Ti Jtm 

and Sb = since f(T,;) — f{Tf). As a consequence, Z[(j)p\ = J dxoF{xo) — 1, 
and therefore ^^[^p] = 0. However, since i7[(/)] generally takes on nonzero 
values outside the subset of physical potentials, its functional derivatives 
at a (/)p can take finite values [30 . As mentioned earlier, those functional 
derivatives are indeed quite remarkable since they are simply related to key 
physical quantities. 



3.4 Functional derivatives of 

3.^.1 Quick comparison with quantum action junctionals 



Action functionals of the form (311, i.e. the exponential of an action, play 
a fundamental role in field theory, which often deals with generating field 
integrals (i.e. traces) of the form |31) . 



e 



#1 £^^^2 (V'l I Po I ^2) r^V^er^ (•^'J 



where is the field, 5Q[tf'] is the action, po is the initial density matrix, 
and (f) is an external source term that couples linearly to ^. In that case, 
the whole exponent in Eq.( |32[ ) is linear in 0, and the successive functional 
derivatives 5^"^/2q[0]/(5(/>(1) . . . 6(j){n) simply insert the field ^ in the field in- 
tegral and generate the averaged (time-ordered) products {54>{l) . . . 54>{n)) of 
the field, alternatively referred to as the time-ordered correlation functions, 
propagators or Green's functions of the theory. In contrast, the classical ac- 



tion functional (31 ) is highly non-linear on the external potential (p since the 
phase-space trajectory x{t) in iS[(/>;a;o] implicitly depends on the external 
potential (j) as well. If, as already mentionned in the previous paragraph, one 
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regards S'[0; xq] as the classical contribution to the total quantum amplitude 
to go from xo to x{Tf), we see that the reduction of the quantum field integral 
from all possible field configurations to the classical path only is responsi- 
ble for the non-linearity in the external perturbation of the classical action 
(31 ). As a consequence, two types of quantities arise when differentiating the 
classical action functional instead of only one, namely the time-ordered 
correlation functions, in quantum field theory. The explicit linear dependence 
oi S on (f> generates the correlations of the phase-space density as in field the- 
ory, while the implicit non-linearities generate terms in volv ing the Poisson 
bracket [•, ■]pb such as the response functions x^'^^j Eq.(21) [35]. 



3.4-2 First and second derivatives. 



In the neighborhood of the external potential can be expanded in 

the series, 

oo „ 

n[(j)o + S(f>] f di...dhn^''\i,...,h)5(t){i)...5(l){h), (33) 

n ^* ^ 



where we introduce the short-hand notation 

f + 00 



dl. 



dTit'{Ti) / dqi / dpi 



V 



Here — i^[0o] and, for n > 1, i?*^") is the n-th functional derivative at 
potential <po, 



r2(")(i,. 



, n) 



= 00 



The latter can be obtained by a systematic perturbation expansion in powers 
of the variations 64> around 0o- We are particularly interested in this work 
in the first two derivatives n = 1,2. For clarity, we just report here the main 
results and give their proofs in appendix [Cj in addition, the differentiability 
properties of are discussed more thoroughly in appendix [d| 

The first functional derivative generates the phase-space distribution func- 
tion / of the dynamics under investigation. 



5n 



= n{r){N{r, p, r)> = i7(T)/(r, p, t(T)) ^ /(l) , (34) 



where we recall that 7T(t) is the window function defined by Eq.(30). Ac- 
cordingly, its second functional derivative. 



(S(2)r2[<?!>] 



'5/(1) 



d4>{l)S^{2) 0=00 (50(2) 



0=00 



x(i,2) 



(35) 
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can be regarded as the linear response function along the closed-time contour 
of the system's dynamics. The direct calculation of x developed in appendix 
[O leads to 



2) - n{T,)n{T^) {nil, 2) - -c{i, 2)) 



(36) 



and consist in the linear combination of two terms. The last term is the 
two-body correlation function 



C{i:2) = {5Af{\)6U{2)), 



(37) 

and, as discussed earlier, comes from the explicit linearity of the action S in 
the external potential. The first term in (36) originates from the non-linerity 
of the action and involves the many-body Poisson bracket [., .]pb, 

7e(i,2) = ^(r4AA(i),AA(2)]^^). (38) 



In Eq. (381, we have introduced the chronologically ordered Poisson bracket 



of two dynamical variables Air) and B{t) as 

Accordingly, x can also be expressed as [33], 

2) = 0{n - T2)x> (i, 2) + 0(t2 - Ti)x< (1,2), (39) 
where for ri > T2 (dropping the window functions for simplicity), 

1 



X>(i,2) = -([AA(i),AA(2)]^^)-^C(i,2), 



and for T2 > ti. 



X< (i, 2) = - 2 ( AA(2)] p^)-'-C{l,~2). 



(40) 



(41) 



3.5 Connection to physical response and correlation functions. 

Remarkably, the function x contains information about both the physical 
correlations and response properties of the physical system, which can easily 
be extracted as follows. 

The physical correlation function C(l, 2) defined in Secj2] is easily ex- 
tracted from X by taking the sum of Eqs.(40) and (41 1 



C(l,2) = -(x>(l,2)+x<(l,2)) , ytut2>to. 



(42) 
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The connection to the physical response can be seen by calculating the vari- 
ation Sf of the phase-space density due to a change in the external potential 
S4>o around 0o- We have, 

(5/(xi,ti) = y"d2x(i,2)5</)o(2) 

ti 



(ix2 / dt2X i^l,ti;X2,t2)S<j)oi^2,t2) 



to 

dX2 



X<(xi, ti; X2, t2)S4)o{x2,t2) 
to 



dx2 / x<(xi,ii;x2,i2)<5(/'o(x2,i2) 



= / dx2 / dt2 [x^(xi,<i;x2,t2) - X^(xi,ti;x2,i2)] '5(?!)o(x2,i2) 

Using Eqs.( |40||4l| ), we find (5/(1) = / dx2 fZ,dt2x'^{l,'2)dM^), with 
X^(l,2) = 0(ti-t2) [x>(l,2)-x<(l,2)] 



^?(ti-t2)^?fe-to)([A^(l),^(2)] 



PS /' 



which is just the retarded response function (21 1 of traditional perturbation 
theory recalled in Sec. |2.2| 



In conclusion, simple linear combinations of the second order derivatives x 
yield to key dynamical properties of the physical system under investigation, 
namely C and x^'^- The simulateneous occurence of both C and two 
quantities of different nature, is remarkable and further discussed in appendix 
[E| with regard to its counterpart in quantum field theory. 



3.6 Generalized BBGKY hierarchy for the extended dynamics 



We consider the statistical dynamics on the closed-time contour governed 
by Hd, defined in Eq.(27|. We assume that at the initial time r^, the initial 



positions and momenta x{Ti) are distributed according to the same distri- 
bution Fo(a:o) as the physical system under investigation (see Secj2j) The 
notations used here for ^1,73,... are those of Sec 2.2 with the real time 
variables ti , f2 • ■ 



simply replaced by ti , r2 . . . . 
The microscopic phase-space density J\f{l) evolves according to. 



1 d 

V(T)dT 



-73(i,2,3)AA(2)A/-(3). 



(43) 



By averaging over the initial state, w e ob tain the evolution equation for the 
phase-space density / defined by Eq. ( 34 1 , 



1 d 

¥{T)dT 



1 



/(I) - {u'"^(l),/(l)} = -73(1,2,3)C(2,3) +Zi(r)/(1), 

(44) 
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in terms of the two-point correlation function C and of 
1 dn{T) 1 



A{t) 



t'{T) dr t'{T) 



5{t - n) - 5{t ~ Tf) = 5{t{T) - to) . (45) 



Following the traditional BBGKY approach, the evolution equation (44) 
would be regarded as the first equation of the hierarchy of evolution equa- 
tions between the successive equal-time cumulants C*-"-*. Our approach allows 
us to instead regard Eq. ( 44 1 as the first equation of a "larger" hierarchy in 
terms of the successive functional derivatives . . . ,n), which combine 

information on both correlation and response functions. The advantage of 
this extended hierarchy is that it can be formally closed at second order 
n — 2 using the procedure discussed below. 

To obtain this extended hierarchy, we note the equal-time property 



X{Xi,t;X2,t) 



-C{Xi,t;X2,t) 



as can simply be seen by setting ti ~ T2 ~ t m Eq. ( 39 1 . The collision integral 
in Eq. ( 46 1 can then be rewritten in terms of x ^nd yields the evolution 
equation 

1 d 



i'(r) dr 



/(i) - /(!)} = -73(1, 2, 3)x(2, 3) + Zi(r)/(1) 



(46) 



which relates the first derivative / = Sf2/S(f) to the second derivative x = 
S'^n/Scf)^ of the generating functional. Successive functional differentiations 
of Eq. ( 46 1 with respect to 4> generate the extended hierarchy between the 
I2(")'s. The first functional derivative yields the evolution equation for x (see 
footnote PS] for details). 



1 d 

¥{T)dT 



Li 



x(i,i')-r'"^(i,2)-x(2,i') 



(47) 



= S2(i, i') + ^73(1, 2, 3)I?(3)(2, 3, i') + A{r)n{r')C{i, V)/zs , 
where we define [36], 

S2(i,i')^{'5(i-i'),/(i)}(i)- 

The mean-field kernel is defined as in Eq.(pO|), 

^™/(i,i')=73(i,2,3)/(2)(5(i'-3). 



(48) 



In the rhs of Eq. ( 47 1 , the second term comes from the direct differentiation 
of the collision integral in Eq. ( 46 1 , and the last term includes the boundary 
conditions. 



23 



Similarly, functional differentiation of Eq.(47) yields the evolution equa- 
tion for J7(3), 



1 



d 



t'ir) dr 



/?(3)(2,3,i) - i:'"^(i,2')r?<^n2',2,3) 



(49) 



= S3(23i) + ^73(i2'3') U2(^n2',3',2,3) + -x(2', 2)x(3', 3) 



IS 



where 



S3(23i) = {5{l - 2),x(i,3)} (i) + {(5(1 - 3),x(i,2)} (1) , 



(50) 



and so on and so forth. The hierarchy thus obtained for the i7*^"^'s suffers 
from the closure problem whereby the equation for Z?*^"-' depends on the next 
order response function ]7("+^). In the following we show that this hierarchy 
can be closed at the n = 2 level. 



3.7 Formal closure of the hierarchy [57] 

3. 1. 1 Effective potential and vertex functions 

The method uses a Legendre transformation of [5Tll38j . a functional 

of the external potentials to obtain a function r[f^] of the phase-space 
densities given by 

(51) 



r[f^] = -n[4>]+ ^di/4i)0(i). 



We shall refer to F as the effective potential] the latter can be regarded as 
the generalization of the Helmholtz free energy used in density functional 
theory [22]. The transformation (51) shifts attention from the potential that 



perturbs the system to the phase-space density that describes its effect. 

The functional derivatives of the Legendre transform r[f^\ at the phase- 
space distribution f,j,~f under investigation, 



r(")(i. 



,n) 



S"r[f^] 



5f^{l)...5f^{fi) 



f4>=f 



are intimately related to those of In statistical field theory, the functions 
_r("^ are usually referred to as vertex functions. For n = 2, we shall write 

r = r(2). 

The first derivative simply equals the external potential, 



5r 



= M^) ■ 



(52) 
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In the special case when 0o = as in equihbrium or relaxation problems, the 
previous relation shows that the physical density / under investigation is an 
extremum of the effective potential F. As such, F plays a role similar to a 
thermodynamic potential in equilibrium statistical mechanics. 



The higher-order derivatives can be obtained by differentiation Eq.(52) 
with respect to (/> and using the chain rule of differentiation. For instance, for 
the second-order derivative, we find 

.2 -J^^ I d2 ''M ^S.a^V) (53) 

6Uil)SUi2) S4>iV) I 5U{l)5U{2) SmHiV) 

where the (5c-function 5c(i — 1') — :jij:^5{l — 1') satisfies 

^(i25e(i-2)/(2) = /(i). (54) 



With short-hand notations, Eq.(53) writes as 



r(i, 2) • x(2, i') = x(i, 2) • r(2, i') = 6,{i - V) , (55) 

and shows that F can be reg arded as the inverse of the response function 



Taking the derivative of ( 55 ) with respect to (f) yields a relation between P^^^ 
and i7(3)^ 

(i, 2, 3) = i')x(2, 2')x(3, 3')/^(^ni', 2', 3') , (56) 

and so on and so forth 1381. 



3. 7. 2 Closure of the hierarchy 



The relation ( 56 1 provides the starting point for closing the extended BBGKY 
hierarchy linking the Firstly, we can now express the collision integral 

(|47| as 



-73(i23)/?(3)(2^3j,) ^ -73(123) [-x(2, 2)x(3, 3)x(i', l')^^'n2, 3, 1') 



= r(i,2).x(2,i'), 

in terms of the memory function kernel 



^(1,1') - -|73(i23)x(2,2)x(3,3)r(3)(2,3,i') 
.!£^3(i23)Mi^ 



(57) 
(58) 



The equation of evolution ( 47 1 becomes 



G(i, 2) - r-/(i, 2) - s{i,~2)\ ■ x(2, i') = S2(i, i') + A{T)n{T')C{i, V)/ts 

(59) 
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where for convenience we have defined 



5(1-2). 



Equation ( 59 1 for the closed-time contour response function x '^iH play here 



a role similar to the Dyson equation in quantum field theory [381126) . In 
the latter case, the equation describes the evolution of the time-dependent 
Green's functions, from which all other dynamical properties including cor- 
relation and response functions can be calculated. Here, in constrast, Eq.(59) 



describes the evolution of the quantity x, defined in Eq.(36), which combines 



in one single quantity, two properties of different nature in classical mechan- 
ics, namely the correlation and response functions (see also discussion in 
appendix 



Secondly, Eq.(59) allows us to derive another relation that expresses F'^^^ 



in terms of S. To this end, we multiply Eq.(59l on the right by the inverse 
— r, which leads to a relation between S and F, 

G(i,i')-i:"-^(i,i') - 

= {Fil i'), /(!)} (i) + zi(r)i7(r')C(i, 2)^(2, i')/^s . 

By taking the functional derivative of the previous equation with respect to 
/ and combining the result with Eq.(57), we obtain the following integro- 



differential equation for the memory function kernel, 

{i;(i,i'),/(i')}(i') 

= |73(i23)s3(23i') 

+ ^73(i23)x(2,2)x(3,3) (^j,ii'2'3') + ^-^j^ 



f4>—f 

+ Zi(r')273(i23)x(2,2)x(3,3)^ [C(i', 2')r(2', 2)] 



(60a) 
(60b) 
(60c) 



The set of equations consisting of the first two equations ( 46 1 and ( 59 ) of the 



hierarchy, together with ( 57 1 and ( 60 ) is closed and is fuUv equivalent to the 
hierarchy of evolution equations for the i?'"-*. Equation (60 1 is a functional 



equation for the memory function that can be used to generate self-consistent 
approximations in terms of 73 and x (see also appendix [f|). 



3.8 The memory function kernel 



Equation ( 60 ) implies that the memory function S can be splitted into three 
parts as 



^(1, 1') = I:^"\X,,X[;t)A{t') + S'{X,,X[;T)SiT - r') + ^^^(1, 1') 



(61) 
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The physical interpretation of the different components of S is discussed in 
detail in the following sections. We just give here a short description of their 
origin. 

The first term 

{i:(°)(Xi,X;;r),/(i')}(i') = ^73(i23)x(2,2)x(3,3)^ [0(1' ,2')x-\2' ,2)] 

describes the effect of initial correlations in the initial state at time to on 
the dynamics of the system at t > tg. This term can be shown to vanish for 
Gaussian initial conditions but, as was recalled in Sec |2.2[ this simplification 
is always invalid here since the phase-space density A/ is always strongly 
non-Gaussian '15|. The neglect of this term has deep consequences related 
to so-called Bogolyubov 's condition of weakening (suppression) of the initial 
correlations and to the transition from reversible to irreversible dynamical 
equations .39; . 

The second term, proportional to the delta function S{t — r'), is singular 
in time and satifies, 

{S^'\X^,X[;T)J{l')}{i') = |73(i23)s3(23i') 

= 73(i23){,5(i'-2),C(i',3)}(i') (62) 

This term is related the instantaneous effect of correlations on the effective 
interactions between particles (at equilibrium, those correlations arc the so- 
called 'static' correlations [22V) This term naturally combines to the mean- 



field contribution Z""^ , Eq. ( 48 ) , such that 



= + ^^^^ - 

renormalizes the instantaneous mean-field effect with instantaneous correla- 
tions effects beyond the mean-field. In the absence of particle correlations 
(Vlasov approximation), S^^") vanishes. 

Finally, the last term describes memory, non-markovian correlation effects 
in the dynamics and satisfies 

{^(-^)(i,i'),/(i')}(i') 

5S{V,2) 



|73(i23)x(2, 2)x(3, 3) (^73(i'23) + ^^^^-^ 



4 Real-time formalism (I) 

The previous approach allowed us to derive a closed, self-consistent set of 
equations to describe the evolution of the single-particle phase-space distri- 
bution function / of the system. However, being expressed in terms of quan- 
tities defined along the closed-time contour, this formulation lacks physical 
transparency and does not appeal to intuition. In the next two sections, we 
shall re-express these exact results in terms of quantities defined on the real- 
time axis, namely the correlation C and response functions x^''^ and a set of 
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memory functions. In the following section, SecjSj we will recast the overall 
approach into a form that leads to the same results directly in terms of phys- 
ical quantities defined on the real-time axis. In the present section, we focus 
on the evolution equations, discuss their physical content, and make contacts 
with previous works. The careful derivation of the evolution equations for C, 
-^R,A (jQjjg ^jjg next subsection is somewhat lenghty; the reader uninter- 



ested in those details may safely skip ahead to Eqs.( 70 ) in Sec 4.2 where the 



results are given and then discussed and compared with previous works. 

4.1 From closed-time to real-time representations 

Throughout the section, we will use the (Pauli) matrices defined as 

'O l\ f -i\ f I 

^.x = I 1 ' = U ' = -1 



4.1.1 Real-time (matrix) representation of x 

The most straightforward transcription of the closed-time contour formalism 
is obtained by mapping the function 1') onto the 2 x 2-matrix 

^^^'^>- [x-+{i,i') x-(i,i') ') ^^-^^ 

defined such that the jfc-component, i,j = ±, is 1') with 1 lying on Ci 
and 1' lying on cw^where c+ (c_) is the upper (lower) part of the closed-time 
contour (see Figjl]). Explicitly, using Eq.([36|, the matrix components are 



X++(l, 1') = I {Tc [7V(1), iV(l')]pB) - ^^(1, 1') 
X- (1, 1') = I {Ta [N{1), N{l')]p^) - JC(1, 1') 
X+-(l, 1') - ([7V(l),iV(l')]ps) - -C{1, 1') 



X-+{lA')^\{[N{l),N{l')]p^) ^-C{IA') 

where Tc and Ta are the chronological and anti-chronological time-ordering 
operators on the real-time axis [io,-l-oo]. This matrix representation has, 
however, the disadvantage that it does not explicitly take account of the 
linear dependence of its components; for instance, 

X++ + X — = X+- +X-+- 



It is possible to chose another representation that not only removes this re- 
dundancy of information but also involves the physically relevant quantities. 
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Indeed we note that simple linear combinations of the Xjk^ yield the corre- 
lation and response functions of interest, 

x" = X++ - X+- = X-+ - X— 
X"^ ^ X++ - X-+ = X+- - X— 
2 ^ 

— C = x+- + X-+ = X++ + X 

IS 

Using these relations, it is straightforward to show that the orthogonal trans- 
formation [40] defined as x = QxQ^ with the orthogonal matrix Q = (1 — 
iay)/\/2 defines the equivalent real-time representation, 



X = QxQ^ =1 , (64) 

which has the remarkable property that it depends on the correlation function 
C and response functions x^'^ only. 



4- 1.2 From closed-time to real-time integrals 



The theory derived in Sec|3] involves integrals over the closed-time contour 
of the form f d2E{l, 2)x(2, 1'). Here we provide useful identities to express 
such quantities in terms of integrals along the real axis (see e.g. |26j.) 

We remark that the quantities x and S defined before belong to the class 
of functions of two closed-time variables of the form 



(65) 



+ f\r)d,{T - r') + I 



r')/>(T,r')+0(r'-r)/<(r,r'), 



where for simplicity we do not write the possible dependence on phase-space 
variables. In matrix for m, / can be represented by / as in Eq.(63) or by 
/ = QfQ^ as in Eq. ( 64 1 , namely 




where the direct calculation of the retarded (R), advanced (A) and correlation 
(C) components in physical time gives the following expressions. 



/«(<, t') = fmt - 1') + e{t - 1') [/> (t, t') - f< it, t')] 

f^{t,t') = fmt - 1') - e{f - 1) [.r{t,t') - f<{t,t')\ 

-fit, t') = 2/W it')5it' - to) + (/> it, t') + /< it, t')) , 



(66a) 
(66b) 

(66c) 



In the special case f — Xi equations (66) give back J^'"^ = x^''^ and — C. 
With f = S, the retarded and advanced terms S^'^ contain the singular 
part that arises from static correlations, while S'~' contains the part 
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due to the initial conditions (see Eqs.( 7l]|72 ) below). Finally, for the 
delta function Sc defined by Eq.(54), we obtain the matrix representation 

4(i,i') = F^<5(i-i>.. 



We are interested in integrals of the kind 



/i(t,t') = * dTf{T,T)g{T,T') , 



(67) 



where / and g arc of the form ( 65 ). It is easy to show that /i(t, t') also belongs 
to the class (651. By splitting the closed-time contour integral into integrals 
over the forward and backward branches, namely §^ dr — J^^ dt — j^^ dt, it 
is straigthforward to derive the following relations, 



h'\t, t') = ■ g"it, t') = / dt"f{t, t")g^{t", t') 

J to 



h^{t,t') = f^-g^it,t')= / dt"f\t,t")g^{t",t') 

JtQ 

h^{t,t') = f-9^{t,t') + f"-g^{t,t'). 
In matrix form, those properties can be summarized as 

Ji(l,l') = /(l,2)aJ(2,l') , Ml,l') = /(l,2)a,g(2,l'), 

with 



f^-g^ 



(68) 



4.2 Equations of motion 

With the help of the previous identities, we can easily re-express the evolution 
equation (59) of x(lj 1') into equations for the physically relevant quantities 
C(l, 1') and ^■'^^'^(l, 1; ) in real time. In matrix form, Eq.(59) writes 



— — Z"-^ G — S"^^ — J ' 



= {^(1-1'),/(1)}CT. 



(69) 



where we recall that the phase-space density /(I) evolves according to Eq.( 18 ) 
In Eq.(69), we introduced the short-hand notations. 



G(l, 2) = -^'5(1 - 2) + {h^{2) + 0o(2),<5(l - 2)} (2) 
Z™^(l,l')=73(l,2,3)/(2)<5(l'-3) 
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Explicitly in terms of the components, we obtain the coupled evolution equa- 
tions for the correlation and response functions [41j . 



d_ 
di 



Li 



dt 

d_ 
dt 



Li 



C{l,l')-J^ dhj dX2Ullt{l,2)C{2,l') (70a) 
= ^ dt2j dX2 2)x^(2, 1') + S{t - to)C{l, 1') 

X^(l,l')-y" dhj dX2Slldh2)x"P,l') (70b) 

= {/(l-l'),/(l)}(l) 

X^(l, 1') - J^^ dh J dX2 Stih 2)x^(2, 1') (70c) 

= {/(l-l'),/(l)}(l) 
Here we have combined the mean- field contribution Z""-'^ to S^'"^, 

i:,Y(i,i') = i:'"^(i,i') + ^'''^(i,i') 

= iJt,6{t - 1') ±e{T{t- 1')) {S>{1, 1') - Z<(1, 1')) 

^ Sl,6{t - t') ±0{T(t- i'))^^(l, 1') (71) 



and, according to Eq.(66c), S'-" is given by 

1') = {Xi,X[ ; t)5{t' - to) + ^ {S> {t, t') + S< {t, t')) 

= {X, , X[ ■ t)5{t' -to) + S^^^''^ [t, t') , (72) 

where in the last expression we introduce notations where is does not appear 
explicitely since indeed isS''^^ and is{I]'' — Z*^ ) a re independent of s (see 
e.g. Eq.(42 ) for a similar property.) With Eqs.( 7l]|72 ), the evolution equations 
(70) can be further details as follows: for all > to, 



dt 



C(l,l') 



(73a) 



dX2Slt{XuX2-M)C{X2MA')- / dt2 / dX2r^(l,2)C(2,l') 



to 



= j dX2i:^^"(Xi,X2;ti)x-^(^i,io;l') + f dt2 J dX2i:^'^^3(l,2)x^(2,l') 

+s{t-to)c{i,r) 

and for all to < t' < t, 

\-Li x^(l,l') (73b) 

J dX2SL{X,,X2;t,)x''{X2,ti,r) - dt2 J dX2S^{l,2)x"i2,l') 
{5(1-1'),/(1)}(1) 
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and for all to < t < t' , 



^ L 



(73c) 



J dX2Stt{Xi,X2;t,)x''{X2,ti,l') + dt2 y"dX2 1:^(1, 2)x^(2,l') 
{J(l-1'),/(1)}(1) 



In the remaining of this section, we discuss general properties of Eq.(70| 
and make contact with previous approaches. The full closure of the set of 
equations for f,x^'^ ^'^d in real-time is completed in Secjsj where we 
transcribe the closure relation (60) in terms of J]^^^ '^ and /, x , C*- 



4.3 Discussion 



Initial conditions. Equations ( |18| and ( 73a||73c ) specify an initial- value prob- 
lem for which t = t' = t^ is the initial time. They require initial conditions 
for /(X,io), x''(^,io;^',io) = -X'^(^,to;^',io) and C{X,to; X't^). Since 
y^RA ixLust comply with Eq.(113l, i.e. 

X^/^(X, io; X', to) - T {5{X - X'), f{X, to)} {X) , 



their initial value is fixed when f{X, to) is given. Therefore, the initial con- 
ditions are completely specified by the given of f{X,to) and C{X,to;X'to). 



form a set of causal equations with 
time history integrals. The presence of those memory integrals is a property 
of the exact statistical dynamics defined in Sec[2] and, despite their formal 
resemblance with (generalized) Langevin equations, they should be distin- 
guished from phenomenological non-equilibrium equations. Thus the com- 
plete information on time reversibility and other symmetries of the micro- 
scopic Hamilton dynamics resides in the time dependent, non-local memory 
kernels: the equations themselves do not single out a direction of time. The 
transition to irreversible dynamics can be induced by approximation, e.g., 
as alluded earlier, when neglecting the dependence on the initial conditions 
(Bolgolyubov's condition). 



Reversibility. The equations ( 73a||73c 



Many-body effects. The original non-linearity of the microscopic equations 
(13) has become separated into several effects [42]. Given a certain level of 
fluctuations, C(l,l'), the term S^^^ can cause them to grow or decay, and 
when viewed as a matrix in its momentum indices, transfer fluctuations from 
of component of C to another. The singular part -|- is related to the 
instantaneous effect of the mean-field and of the correlations, respectively, 
while the regular part describes non-Markovian (i.e. delayed) correlation 
effects. Another effect contained is the renormalization of the spectrum 
of fluctuations by the many-body interactions. 
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The term on the left-hand side of Eq. ( 73a I can be interpreted as follows 
(a more complete discssion is given in Sec 4.6 j Given a fluctuation 5N of the 
density around its average /, let S(j)(l) be the potential defined such as 

In other words, we imagine that SN can be created by disturbing the system 
by the external potential Scf). From Eq.([70|, we have 



(G - - S^) (1, 2) • SN{2) = {0(1), /(I)} . 

Multiplying the previous equation by 6N{1') = x^(l', 3)<5</)(3) = x^(3, l')5(j){3) 
and ensemble averaging yields 

(G - S^^f - S"") (1, 2) . G(l, 1') = mm3))J{l)} ■ x^(3, 1') . 

Comparing this with Eq.(70), we may interpret as an effective random 



source of potential fluctuations. Since, as we shall find later, includes 
terms which are quadratic in G, this term can be interpreted as a source 
of non-linear noise which drives the fluctuations G, and couples modes of 
different wave numbers. 



More on those properties will be discussed in Sec 4.6 In the next two sub- 
sections, we compare the results derived so far with previous works, starting 
with fluids at equilibrium and then with non-equilibrium kinetic theories. The 
comparison with those previous works is not thorough and would certainly 
deserve a more detailed analysis. 



4.4 Equilibrium limit, fluctuation-dissipation theorem and detailed balance 

When considering a system at equilibrium, we generally implicitly set the ini- 
tial time to = — oo. As a consequence, statistical averages like the correlation 
and response functions are invariant under time translation, e.g. |43j 

C{h,t[) ^C{h + s,t[ + s) , VseM, 

and are functions of the time difference t = ti — t'l, i.e. 

C{l,l') = C,giX,,X[;t,~t[) , yt,,t[. 

As recalled in Sec. [2j a fundamental property of equlibrium fluids is the 
fluctuation-dissipation theorem, which simply relates C^q and Xeq^ 

X^qiX,X';t)~Xe,{X,X';t) = ~p§-CeqiX,X';t) . (74) 



with the inverse temperature j3 — 1/KbT. By rewriting the evolution equa- 
tions (73) in terms oit = ti — t'-^, it is straigthforward to show that these equa- 
tions are consistent with the fluctuation-dissipation theorem if and only if 
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the equilibrium memory function kernels satisfy a relation similar to Eq.(75) 
[46], 

E^^iX, X'; t) - E^^{X, X'- 1) = X'- 1) . (75) 

This relation describes detailed balance of the collision mechanisms that con- 
trol the occurrence of spontaneous fluctuations and their damping. 

It turns out that the evolution equations obtained by setting Iq — —oo and 



t = ti = t'l in Eqs.(73) are rather cumbersome to manipulate. In the study 
of the dynamics of equilibrium correlation functions, it is more customary to 
rather work with the following explicit expression of the correlation function, 

Cm{X, X'; = C{X, t; X',0) = {6N{X, t)SN{X', 0)) . (76) 

This expression is one possible representation of Ceq{X, X' , t) (corresponding 
to s = in C(t + s, s) = Ceq{t).)\ contrary to Ceq, Cm favors an initial time, 
io = = 0. Previous works on the equilibrium kinetic theory of fluids and 
particular Mazenko's theory [12] discussed in the introduction are theories 
for CM{t)- In our approach, the evolution equation for Cm can simply be 



obtained by setting both the initial time to and the time t' to zero in Eq. ( 73a I . 
We immediately obtain Vt > 0, 

^-Li(X) CM{X,X';t) - J dX2Stt{X;Xr,t)CM{X2,X',t) 

~ lo '^^ I ^^2^'^(^,i;^2,t)CM(^2,^';i) = 0. 

This evolution equation has the well-known form of memory function equa- 
tions in the equilibrium theories [43] . In particular, it corresponds to Mazenko's 
memory function (j) defined in |12j with the identification, 

El,{X,X-t) = -cl>\X,X) 
S^{X,t;X',t') = -(l)''{X,X';t^t') . 



4.5 Contact with other non-equilibrium kinetic theories 

In this section, we explore the correspondence with previous non-equilibrium 
kinetic theories. First we make contact with some of the most popular closures 
of the BBGKY hierarchy. We shall see that those closures systematically 
neglect the effects of the memory function kernels S'^''^ . Then we compare 
our results to the theory of Rose |15l mentioned in the introduction. 
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Closure 


A 


B 


C 


D 


E 


<5 








Landau 


X 




















LBG 


X 


X 

















Boltzmann 


X 





X 





Eq.( 


79 










Book-Frieman [1?] 


X 


X 


X 





Eq.( 


79 











Tabl e 2 Left side: terms kept (x) and dropped (0) in the second BBGKY equa- 
tion (77 1 for the Landau, Boltzmann, Lenard-Balescu-Guernsey and Book-Frieman 
closures. The latter, which was proposed to improve the close collisions in the 
Lenard-Balescu-Guernsey equation, simply amounts to setting the three-body cor- 
relation (73 equal to zero. The same closures are expressed in terms of the cumulants 
C'-"-' in tablehl Right side: in the approach developed in this paper, the approxi- 
mations for Cand D correspond to the approximations for the memory functions 
j^c5,A,c gjygjj jjj ^jjg jg^g^ three columns. 



4-.5.1 Reduction to popular BBGKY closures 

The BBGKY hierarchy is a hierarchy for the equal-time correlation functions 
gn (for convenience, we give a short reminder of the BBGKY hierarchy in 
appendix|Aj) Traditional closures of the BBGKY hierarchy are performed at 
the level of the evolution equation for §2, which reads 



52(l,l') = iii'/(l)/(l') (77) 



- Jd2L,2f{l)92{l', 2) + Li.2/(1').92(1, 2) + [^12 + Li'2] /(2).g2(l, 1') 



+ Lii'92 



(!,!') + y"d2(Li2 + Li,2) 53(1,1', 2), 



where, throughout this sub-section, we write 1 = (r, p) and all quantities 
are evaluated at the same time t = ii They rely, among other 

hypothesis, on an assumed ordering of the correlations .g„ in terms of an 
adequately chosen small parameter A 

gn = o(A") . 

The Landau, Boltzmann and Lenard-Balescu-Guernsey closures give evolu- 
tion equations for / that are second order in A. The Landau closure is valid 
for so-called weakly coupled systems and A is the dimensionless strength of 
the potential (w(r) — \v{r) in Eq.Q.) The perturbation parameter in the 
Boltzmann equation is the density parameter A — nl^ where n is the par- 
ticle density and Ic is the correlation length. The Lenard-Balescu-Guernsey 
equation is an equation for weakly-coupled plasmas and A is the so-called 
plasma parameter A = q^/aksT {q is the particles charge, a = n~^/^ is the 
interparticle spacing.) In the Lenard-Balescu-Guernsey equation, the plasma 
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is represented as a weakly coupled system in which the collisions are due to 
the interactions via an effective, dynamically screened potential. For those 
popular closures, table [2] shows the terms that are kept (x) and dropped (0) 
in the evolution equation ( |77[ ). 

In our approach, the equation for 172(17 1') = 1') ^ <^(1 ^ can 
be easily obtained from Eq.(73a) as 

d 



dt 



[Ll+Ly) 



.g2(l,l')-iii'/(l)/(l') 



(78a) 



- Jd2L,2f {1)92(1', 2) + Lr2/(1')32(1, 2) + [L12 + Li'2] /(2)g2(l, 1') (78b) 



^2^:^(12)32(2, 1') + U\V2)g2{2, 1) (78c) 
d2^^(12)52(2, 1') + 52(1, 2)^^(21') (78d) 
+ J d2S^{12)x^{2, 1') + x^(l, 2)^^(21') , (78e) 
where t — t' — t2 is implied. We note that in the rhs, the terms on lines 



(78a I and (78b I, which are just the A and B terms of Eq.(77), come from the 
mean- field term , while all th e ot her terms come from memory functions 
jjR,A,c' ^ The term C + D of Eq.(77) are replaced by the last three lines in 



Eq.(78l 



In order to appreciate better the physics contained in those memory func- 
tion terms, we translate the tradiational closures discussed before in terms 
of approximations of Eq.(78); the results are summarized in table [2] 

Both the Landau and Lenard-Balescu-Guernsey closures discard all the 
terms involving the memory functions, which amounts to setting JJ^''^''^ = 

only and Landau-Balescu-Guernsey discards the 
.) This simply tells us that memory functions will 



(Landau retains Eq.(78a 



terms (78c) through (78e 



account for effects not accounted in those closures, for instance by adding 
the effects of static (S^) and dynamic {S'^) correlations on the effective, dy- 
namically screened interparticle potential through which particles mutually 
interact. 

The Bol tzma nn and Book-Frieman equations are recovered by (i) dis- 
carding Eq.(78b|, (ii) setting ~ S'^ — 0, and (iii) using the following 



approximation for S , 

S^ill') = 73(123)(5(1' - 2),5(1' - 3) = 



aw(l-l') d&{l-l') 



dr 



dp 



(79) 



which, together with Eq.( 78c I yields the term C of the second BBGKY equa- 
tion (ITTI) that describes 



Dare two-particle interactions (uncorrelated binary 
collisions) . The approximation ( 79 ) can be obtained using Eq. ( 62 ) together 
with the approximation C(ll') w S{1 — l')/(l) (see table [T]) Indeed, we 
obtain 



{i:^(i,i'),/(i')}(i') 



73(123){5(1' 
73(123){5(1' 



3),<5(l'-2)/(l')}(l') 
3)5(l'-2),/(l')}(l') 
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which imphes Eq.( |79[ ). 

In summary, this quick comparison reveals that inclusion of any coher- 
ent approximation for the memory functions E^'^^^ will bring in physics 
information that is systematically neglected in traditional closures. 



4-5.2 Connection with MSR and Rose's theories 



General remarks. The present theory retains two main imprints of quantum 
mechanics, namely the classical amplitude ei"^ with the action parameter s 
and the closed-time contour idea of Schwinger. In contrast, MSR [14, is based 
on the concepts of canonical quantization using field doubling via conjugate 
response field and, in Rose |15l, on the second quantization using the occupa- 
tion number representation in phase-space. Thus, MSR and Rose's theories 
recast the classical problem in a form identical to quantum mechanics, while 
our approach should be seen as a classical limit of the quantum theory. 

Another difference is in the nature of the external coupling . In MSR 
and Rose, as usually is the case in quantum field theory |31j, the coupling 
corresponds to a source term in the evolution equation of the field, while here 
the coupling is an additional term in the Hamiltonian. The response function 
in MSR and Rose thus describes the response to an infinitesimal source of 
particles while in the present theory it describes the response to a disturbance 
produced by an externally applied force field. Since all three theories use 
similar renormalization techniques based on a generating functional and a 
Legendre transform, the equations obtained look alike. As we shall see, our 
results encompass the results of Rose's renormalized kinetic theory |15) for 
any initial conditions. The ingenious formal devices introduced by MSR and 
Rose naturally arise in the present formalism. 



Technical comparison. Briefly speaking, Rose's theory is a theory for the 
correlation function C and the response function R that measures the linear 
response of the phase-space density to an infinitesimal source term rj in the 
field equation (18), i.e. 



d_ 



Li 



1 



/(I) - {u'-f{l)J{l)} - -73(1, 2, 3)C(2, 3) = 77(1) (80) 



and 



i?(l,l') 



Mill 



The function R can also be regarded as the Green's function of the linearized 
version of Eq.( 18 ). 

On the other hand, in our theory describes the linear response to 
a variation 0o + ^(f'o of the external potential. In the presence of S(j)o, the 
evolution equation (18) can be regarded as Eq.(80) with the source term 



r;(l) = Wo(l),/(l)} 
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which has the particularity of conserving the number of particles (since the 
momentum integral vanishes.) In the linear regime, we can write 

5/(l)-x''(l,2)<50o(2) 

= i?(l, 2)7j{2) = -{R{1, 2), /(2)}(2)<50o(2) 

and therefore, 



x^(i,i')--{i?(i,i'),/(i')}(i') 



(81) 



When introduced into the evolution equations ( 70 ) for t^^'^ and C, we find 



dt 



i?(l, 1') -Jd2 (1, 2)i?(2, 1') = ^(1 - 1') 



(82) 



and 



C(l,l') -yd2i:«,(l,2)C(2,l') (83) 
d2 2), /(2)}(2)i?(l', 2) + d{t - to)C{l, 1') 



These equations are equivalent to the (Dyson) equations (63) and (64) of [TS] 
derived by H. Rose, with the correspondence 

^_+(i,2) = i:,« (1,2) 

i;__(l,2) = {i;^^-«(l,2),/(2)}(2) (84) 



= S' 



CO 



The last equation is the consequence of the fact that Rose assumes Gaussian 
initial conditions. Although as noted by Rose in his appendix A, formalisms 
developed by Deker [JS] and more recently by Andersen [IS] exist to extend 
their validity to general initial conditions, we have not performed the exten- 
sion ourself yet. The closure relations, i.e. the functional equations for the 
memory functions, in Rose and our formalism (see Sec[5]) look quite different 
(simply because they involve different starting points and ingredients) and 
are thus difficult to compare. Moreover, Rose points out that the closure 
relations (52)- (62) of his paper are "exact if and only if the initial condi- 
tions satisfy a Wick-type theorem" , i.e. are Gaussian. Nevertheless, assuming 
Gaussian initial correlations in our formalism, the lowest order approxima- 
tions in both approaches are the same (see Secjsj) At this point, it is not 
clear which approach is more handy when it comes to practical calculations. 



38 



4.6 Formal solution for C(l, 1'). Dependence on initial conditions. 



In this paragraph, we show that it is possible to solve at least formally the 
evolution equation Eq.( |70a[) for the correlation function C(l, 1') in terms of 
the response functions and Green's function R, respectively. When in- 

troduced in the colhsion integral Eq.(18), this solution may be used to build 



a kinetic equation for the distribution function /. In its simplest approxima- 
tion, this approach corresponds to the quasi-linear theory used for instance 
in plasma physics to derive the Lenard-Balescu equation i5Q\ . The solution 
given below is exact and takes care of the initial conditions. 



General expression. Using the definitions of Sec 4.1 it is straightforward to 
show that the inverse F of the closed-time response function x defined by 
Eq.(55l can be represented as [5T] . 



r = 





pR 2_pC 



where the components and merely are the inverse of the retarded 
(advanced) response function amd x"^, respectively, i.e. 



r«(i, 2) . x"i2, 1') = X^'ih 2) • r«(2, 1') = S{1 - 1') 
r-4(i, 2) . x^i2, 1') = x^(i, 2) • r-4(2, 1') - S{1 - 1') 

and F'-' is related to the correlation function as, 

r^(l, 1') = -r^(l, 2)C{2, 3)r^(3, 1'). 



(85) 



(86) 



(70a 



Equation ( 86 1 implies the following formal solution of its evolution equation 



C(l, 1') = -x^'(l, 2)F^ (2, 3)x^(3, 1') = -x^(l, 2)x''(l', 3)r^ (2, 3) (87) 

in terms of x^'^ an d of the quantity F'~^ . Another expression can be obtai ned 
using the relation (81) between x^ the Green's function R in Eq.(87) 
and yields 



C(l,l') = i?(l,2)i?(l',3)5(2,3), 
where we introduce the symmetric kernel, 

5(1,1') = -{{^''(1,1'),/(1)},/(1')}- 



(88) 



(89) 



Equation ( 86 1 can easily be verified by substitution into the evolution equa- 



tion (70a), while Eq.(88) is solution of 



d_ 
Ft 



C{1, 1') - 2 dh j dX, 2)C{2, 1') (90) 

dX2 5(1, 2)i?(l', 2) + 5{t - to)C{l, 1') , 



dU 



to 
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obtained by rewritting the rhs of (70a) with Eq.(81) IH^]. They may be con- 
sidered as a generahzed form of the fluctuation-dissipation theorem in the 
sense that they express the overaU level of density flucutations C as balance 
between forcing (through I^*^ or S) and dissipation (encapsulated in ^^'^^ or 
R) [13] ■ In the next paragraph we show how Eqs.( 86p8 ) can be further split 
into more fundamental components. 



Detailed expression. Introducing Eq.(86l into Eq.(70a) leads to the following 
relation between F'-^ and , 

{r^{i, 1'), /(i)}(i') - ^^(1, 1') + Sit - to)C{i, 2)r^{2, 1') , (91) 

which, as proved in details in appendix [Gj imply the following structure for 

r^(i,i'), 

r^{i, 1') = r-9(i, 1') + Sit' - to)r'{x^,x[;t) + s{t ~ to)r'{x[,x,;t') 

(92) 



+Sit-to)6it' -to)r'-'>\x,,x[) 



where r^<^9-Sfi g^j,g given by 



and 



{r-«(i, 1'), /(i)}(i') = I {s>{i, 1') + s<ii, 1')) 

{r\X,,X[;t), /(1)}(1') = tsS('^\X,,X[;t) 

{ {r(°) {x^,x[)j{x^,t,)}, f{x[ ,to)} = Co{x^,x[), 



(93a) 
(93b) 



where Co{Xi, X[) = C{Xi,to; X[,to) is the iniital value of the correlation 
function. Similarly, equation ( 89 ) implies 



5(1, 1') = 5''^s(l, 1') + S{t' - to)S\Xi,X[:t) + S{t - to)S\X[,X,;t') 

+5{t-to)Sit' -to)S\X,,X[) 



with 



S^^^il, 1') = f {^>(1, 1') + S<{1, 1'),/(1')}(1') 
S\X,,X[;t)^is{s^"\X,,X[;t),f{l')]{r) 
S°{X^,X[)^Co{X,,X[) 



Substituting those expressions into Eqs.( 86][88 ), we obtain expressions that 
describe the temporal evolution of the correlation function from its initial 
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value C{Xi,to;X[,to) to C{Xi,t; X[,t') at times t,t' > to, 



C(l,l') = // dX2dX3x''(l;^2,to)r(°'(^2,X3)x^(X3,to;l') (94a) 



(94b) 



dX2dX3x"il,2)r\X2,X3-t2)x^{X3,to;l') 
+ [ dhff dX2dX3x''{l;X2,to)r\X2,X3;t3)x''{3;l') 
+ fdt2f dh [[ dX2dX3x''{l,2)r-'%23)x^{3,l') 

Jto Jto J J 

=jj dX2dX3R{l;X2,to)R{l';X3,to)CoiX2,X3) 
+ ^ dt2 1 J dX2dX3R{l, 2)i?(l'; Xa, to)S\X2, X3; ia) 
+ ^ dia^y dX2dX3R{l;X2,to)R{l',3)S\X2,X-i;h) 
+ fdt2f dh jj dX2dX3R{l,2)R{l' ,3)S'-^s^23) . 

We can shed light on the different components of Eq.(94) by considering 
the dynamics of the density fluctuations Z\iV(l). An equation of motion for 
Z\iV(l) = iV(l) — /(I) is easily obtained by substractingEq.( 18 1 fromEq.(16). 
By substracting the quantity S^{1,2)6N{2) on both sides of the resulting 
equation, we find 



- - Li^ SN{1) - £'"■'■(1, 2)SN{2) - 2)6N{2) 



- ^73(1, 2, 3) {N{2)N{3) - C(2, 3)) - r«(l, 2)^7V(2) = SJ^{1) , 

which can formaly be integrated in terms of the Green's function i?(l,l') 
(recall Eq.Q) as 



dN{l)^ J dX2R{l;X2,to)SN{X2,to)+ J dt2 J dX2R{l,2)5F{2) (95) 

Thus the microscopic phase-space density fluctuation is written as the sum of 
two terms. The first term represents the propagation in the fluid of the density 
fluctuation that would be caused by a small initial perturbation 5N{X2,tQ) 
in /. Roughly speaking, it describes how a microscopic flucutation propa- 
gates in the fluid on average and neglecting non-linear effects. The second 
term in Eq. ( 95 1 describes the corrections to this interpretation that manifest 
in the equation of motion in the source term 5F and in the solution ( 95 ) 
through non-local effects in both space and time. Using Eq.(95l to build the 
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correlation function C(l, 1') — {6N{1)5N{1')) lead to the four contributions 



in Eq.(94b) with 



S\Xi,X[;t) - {STil)5N{Xlto)) (96) 



Thus, the first term in Eq.(94b) describes the contribution obtained when 
propagating the initial fluctuations as if they were small, independent pertur- 
bations on the fluid. The three other terms represent the corrections to this 
picture, which allows us to express the components of the spectral function 
S in terms of correlation functions between the noise term SJ^ and the initial 
fluctuations 6N: S'^'^^{1, V) is the autocorrelation function of the noise term 
ST while S^'^^ describes the influence the initial conditions. 



Equation ( 94b I exhibits the influence of initial correlations in the ini- 
tial state; the latter can be important when considering ultrafast relaxation 
processes when t, t' approach tg- The influence of initial correlations is often 
neglected in kinetic theory by invoking the condition of complete suppression 
of initial correlations (a.k.a Bogolyubov's condition [35]) according to which 
all initial corraletions existing at tg are damped at a sufflciently large time 



t — » tcor- Inspection of Eq.(94b| sugges ts t hat tcor is determined by 
the decorrelation time of {Xi , X[; t) in Eq. ( 96 1 and by the properties of 



1'). Then for times t — to >> tcor, we set S^^= and 

Co{X,,X[) - Ct''\X,,X[) = /o(r,p)<5(r - r')<5(p - p') 

where /o(i", p) = ( J2f=i (r — S {p — Pj)) is the initial single-particle 
distribution function such that 

C(l,l') ^JJ dX^dX3R{l;X2,to)R{l';X3,to)Ct''\X^,X3) 
+ j dt^j dt3 If dX2dX3R{l,2)R{r ,3)S'-'3{23) . 

J to J to J J 

When dynamical correlations are alltogether discarded {S^'"^''^=0), then 



C{1,1') dX2dX3R"'Hl;X2,to)R"'f{l';X3,to)Ct'''{X2 



with - Li(l)] 1') - 2)R{2, 1') 5{1 - 1') corresponding to 

the linearized Vlasov equation, which in plasma physics is known as the 
quasi-linear theory and leads to the Lenard-Balescu-Guernsey equation. 



5 Real-time formalism (II) 

The previous section focused on the derivation of the coupled evolution equa- 
tions for the correlation and response functions in terms of the JJ^'^^'^' Here 



we concentrate on translating to real time the closure relations (57) and 
(60) to obtain the equivalent relations for the real-time components Z'^'"^'"-^. 
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Although it is possible to directly translate Eq. ( 58 ) in terms of the matrix 
components X'^^'^''^, XJ^'^''~' and of the components of the tensor obtained 
from the different components of the three-point vertex F^^^^ [551151] . the 
manipulations are not trivial and we prefer the approach described below. 
A considerable advantage of this alternative approach is that it allows us 
to introduce an equivalent formulation of the closed-time contour approach, 
which directly generates the correlation C and response functions x^'^- 



5.1 Real-time representation of the generating functional 



Given a potential 0(t) on the closed-time contour, we define the physical 
component 4'p{t) and non-physical (pAit) on the real-time axis as 

1 



0p(X,i) = -[0(X,T+)+(/)(X,T_)] 
jA{X,t)^C^{X,T+)^C^[X,T^) 



(97) 



where t = t(T_|_) ~ t{T_). In this representation, physical potentials are simply 
characterized by 4>a = 0. Similarly, we define the 'physical' and 'non-physical' 
components of the phase-space density as [55] , 



U.p{X,t) = \ [AA(X,r+) +AA(X,r_)] 
^ A/-^ (X, t ) = AA(X, r+ ) - AA(X, T_ ) 



(98) 



For a physical potential, the phase-space density is equal on both sides of 
the contour and therefore Ma (A, i) = and Mp (A, t) equals the phase-space 
density of the system under investigation in the external potential (pp. With 
these definitions, the coupling term in the action (291 writes 



to 



dAi0p(l)AA4(l)+04(l)AAp(l), 



(99) 



which shows that (j)p linearly couples Ma while pA couples to Ma', accord- 
ingly, variations with respect to (pA will generate the physical phase-space 
distribution under investigation. 

Equation ( 97 ) defines a one-to-one mapping between the potentials (j) and 
the vector of potentials {4>p, 4>a)- As a consequence, the generating functional 
can be replaced by the action functional Pa] defined as 



(100) 



In terms of the new variables, the Taylor expansion (33) around the external 
potential 0o becomes. 



^1. 



,12 



] = E fdl.--djfdl'...dk' (101) 
') SPp{l)...S<j>p{j)PA{l')...5<l>A{k') 



j,fe=0 • 
.2(1. ..n;l'. 
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where 

^2i...i2...2(i---j;i'---fc') 



%(l)...50p(j)50A(l')---'50zi(fc') 



1=0 



Explicit expressions for these derivatives can be obtained either from the 
J?*-"-* [0] or from a systematic perturbation expansion as shown in appendix 

The first derivatives generate the phase-space density /(r,p,f) under in- 
vestigation, 



5Q 



%(1) 



= 



5n 



Po,<pA=0 



= /(!)■ 



bp=4>o,'t>A=0 



Most remarkably, the second (partial) derivatives directly generate the re- 
tarded and advanced response functions and correlation functions, 



{ r2ii(l,2) = 
^2i2(l,2) = x-^(l,2) 
r22i(l,2)-x''(l,2) 

1222(1,2) = -C(l,2) 



(102) 



We therefore obtain the noteworthy result that a generating functional de- 
fined in terms of real-time quantities exists that generates the phase-space 
distribution, the correlation function C and response functions x^'^- 

We now continue the clos ure p rocedure by adapting the steps of SecjS] 
to the new action functional ( 100 1 . We define the effective potential as the 
Legendre transform, 



Its deriatives at fp — f, Ja = define the vertex functions 



(103) 



r 



1 . . 

j til 



_12. . 



.2(l...j;l'...fc') = 



Sfpil)...SUU)5fAiV)...SfA{k') 
The first derivative generates the external potential ^o, 



6r 







sr 



fp = fjA=Q 



SfAil) 



0o(l) 



fp = fjA=Q 



while, as shown in the appendix H.l), the second-order derivatives directly 
yield the r^.-^.c^ j_g_ ^ 



r Ai(l,l') = 

ri2(i,i')-r^(i,i') 
r2i(i,i') = /^«(i,i') 

1^22(1,1') --^^(1,1') 



(104) 



IS 
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Higher order vertex functions can be obtained by taking the successive deriva- 
tives with respect to the potentials with the chain rule of differentiation. The 
direct calculations become rapidly cumbersome with the growing number of 
indices and variables and, in appendix [Hj we present a graphical method 
that allows their calculation in a much more economical manner. We just 
give here the results for those derivatives of importance in the present work, 
namely 



X?22i(123) = I222(ll)^22i(22)l2i2(33)ri2i(123) 
+ I22i(ll)^222(22)l2i2(33)r2ii(123) 
+ I22i(ll)^22i(22)l2i2(33)r22i(123) 

-r2222(123) = [l22i(ll)^22(22)l222(33)r2n(i23) 
^22i(ll)^22i(22)l222(33)r22i(123) 

^2l(ll)^2l(22)l22l(33)r222(123) 



c.p 
c.p 



(105a) 



(105b) 



where c.p. stands for cyclic permutation of the triplet (11,22,33). These 
equations are the transcription in terms of real-time quantities of the equation 
( 56 1 derived before. 



5.2 Real-time representation of the closure 

The same steps described in Secs |3.6p?f| to derive the closure can be followed 
with the evolution equations for N^a as a starting point. The main steps are 
given in appendix |l] and can be summarized as follows. From the evolution 
equation for Mp^A, we first obtain the evolution equation for fn A by aver- 



aging and for x^'^ and C by functional differentiation using Eqs.( |102[ ). The 
collision operators in those equations involve the third order derivative s fig k 



and are rexpressed in in terms of memory function S^'^''^' usi ng E qs.( 105 ). 
The overall procedure directly leads to the evolution equations (70 1 for x"'^ 
and C together with the following closure relations 



^''(1,1')= o73(123) 



(5C(2,3) 



r-4(l,l') = zs73(123) 
^^(1,1') = ?73(123) 



'5x^(2,3) 



S.fA{l') 

<5C(2,3) 



fp-jA=0 

fp;fA=0 



(106a) 
(106b) 
(106c) 
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or, after performing the functional differentiations, 



^^(1,1') 



1:^(1,1' 



--73(123) [X221 (22) r22i (33)^221 (2, 3, l') (107a) 
+2/222(22)/?2i(33)ri2i(2,3,l') 

-is73(123)bi2(22)r22i(33)ri22(2,3,l') (107b) 
+r2i2(22)X?22(33)Ai2(2,3,l') 

-^73(123) [/?22(22)/222(33)Ai2(2, 3,1') (107c) 
+21221(22)1222(33)^212(2, 3, 1') + I22i(22)l22i(33)r222(2, 3, l') 



These equations for S 



R,A-C 



are the transcription in real-time of the closed- 



time contour formulas (58) and (57). The set of equations consisting of the 
evolution equations (18 1 and (70) together with Eqs.(106l is formally closed. 

As in Sec 3.7 the evolution equations (70) can be used to derive addi- 
tional relations between the memory and vertex functions to derive functional 
equations for the memory functions £'^^^^c. The details of the calculations 
are given in appendix [ij The resulting closure relations can be expressed in 
a variety of equivalent forms in terms of the quantities C, x^'^^ jjRA,c 
and S encountered before, and we choose to give them in the following form 
(which is most directly comparable with previous works, e.g. |15|): 



^^(1,1') = 73(123)73(231')i?(22)C(33) -)-73(123)i?(22)C(33) 

+ i73(123)i?(22)i?(33)^ 



(108) 



and (see Sec 4.6), 



5-^1,1') = ^73(123)73(1'23)C(22)C(33) + ^j,{123)C{22)Cm^-^^^ {109) 
+ zs73(123)x^(22)C(33)^f^ +z.73(123)x^(22)x^(33)-'^^''^^'^) 



SfA{2) 



and 



S\X,,X[;t) - -73(123)x«(22)x^(33) 



IS 



5/^(3) 



c{x[,to-A}r^{i2) 



In the lowest order approximation, 

S^'il, 1') - 73(123)73(231')i?(22)C(33) 
5'^^s(l,l') = i73(123)73(l'23)C(22)C(33) 
= 
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we recover the famous Direct Interaction Approximation (DIA) derived by 
Kraichnan [58] in the context of fluid turbulence and later extended to the 
Klimontovich equations by Dubois and Espedal [5U1IT5] . Further discussions 
on the application of those closure relations to physical problems would cer- 
tainly require substantial additional work, which is beyond the scope of the 
present paper. 



6 Summary 

We have presented the foundations of a theory to derive closed, self-consistent 
approximations for calculating the statistical dynamics of classical Hamil- 
tonian systems, which can describe both equilibrium and non-equilibrium 
states. The theory, which unifies and encompasses previous results for clas- 
sical Hamiltonian systems with any initial conditions, can be regarded (i) as 
the classical mechanical counterpart of the theory of non-equilibrium Green's 
functions in quantum field theory and (ii) as a generalization to dynamical 
problems of the density functional theory of fluids in equilibrium. The present 
approach avoids many of the complications inherent to previous works |49j ; in 
particular the theory is valid for any initial state and treats equilibrium and 
non-equilibrium states in a unified manner. Interestingly, several of the key 
ideas and tricks of the previous theories (e.g, the need for both correlation 
and response functions, operator doubling, imposed commutation relations 
and other causality constraints...) arise naturally here albeit in a different 
form (e.g., the operator doubling in MSR is replaced by the coupling to 
physical and non-physical phase-space densities Mp^/\, which generate both 
response and corrrelation functions..) 

The main ingredients of the present approach are the following. Given a 
Hamiltonian system in an external potential 0o £^nd in the statistical state 



described by Fq at initial time toi the effective action functional ( 100 ) defined 

as 

^[4'pAa\ = -sln2^[0p,(/>^] , 

where s is a fixed parameter with the dimension of an action, contains all 
information about the dynamical properties of the system at times t > to- 
In particular, its first derivative at (j}p — (pQ, (p/^ — gives the single-particle 
phase-space distribution function, 

with 1 = (r, p, t), and its second derivatives generate the correlation function 
and retared and advanced response functions, 

x-(M) 

X^(l,2) 
C(l,2)/zs. 



J(Ap(l)(5</.4(2) 
J04(l)<504(2) 
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Higher-order derivatives systematically generate combinations of products 
and Poisson brackets of the fundamental field N{1). Using traditional renor- 
malization techniques involving the Legendre transform of the action function 
^[(f'pT 4'a\, a closed, self-consistent set of equations of motion is derived for 
the single-particle phase-space distribution function /, the correlation func- 
tion C = {SfSf), the retarded and advanced density response functions x^'"^ 
to external potentials, and the associated memory functions JJ^'"^''^ , 



d_ 
dti 



/(1)-{«"^(1),/(1)}=-73(1,2,3)C(2,3). 



and 

di 

d_ 
di 



Li X«(l, 1') - d2 2)x«(2, 1') 



^ L 



d2^« (1,2)C(2,1') 



with the memory function kernels satisfying 



d2i;^(l,2)x^(2,l') 



<5x«(2,3) 



IS 



r^(l,l') = ^73(123) 



SfA{V) 
(5C(2,3) 



'5/^(1') 



fp;fA=0 

fp-jA=0 
fp-jA=0 ' 



(110) 



Moreover, the memory functions satisfy functional differential equations, e.g. 
(108 1, that may be used systematic approximations. 

In conclusion, the purpose of the present paper was to lay down the 
foundations of a theory for the construction of renormalizcd kinetic equations 
of classical systems of particles in and out of equilibrium. We hope that these 
results will serve to (re)stimulate further research on this challenging topic. 
Of course further works calls for applications of the approach to realistic 
problems; work along those lines is under way. 
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7. The search for the generalization of the Boltzmann equation has a long and 
rich history [^. It is only after about 1946, under the impetus of Bogoliubov's 
ideas, that concentrated efforts to generalize the Boltzmann equation systemat- 
ically to higher densities was undertaken. Starting with the Liouville equation, 
a method was worked out, notably by Bogoliubov, Choh and Uhlenbeck, Green 
and Cohen and others, which amounts to a systematic inclusion of many-body 
collisions through a density expansion similar to the virial expansion for comput- 
ing equilibrium properties of dense gases. In the mid 60's, it was realized that the 
density expansion of the collision operator was plagued by divergences. While it 
may in principle be possible to deal with divergences by isolating and resumming 
the most secular terms of the virile expansion, the exceeding complexity of the 
procedure has severely restricted the usefulness of the density expansion in un- 
derstanding the properties of even moderately dense gases - let alone dense gases 
and liquids. 

8. The appearance of such divergences has several remarkable consequences. 
Firstly, they imply that the density dependence of the transport coefficients is 
non-analytic (e.g. it has logarithmic contributions). Secondly, they suggest that 
the relation between applied gradients and induced hydrodynamic fluxes in non- 
equilibrium thermodynamics is non-analytic in character, which leads one to won- 
der on the proper macroscopic description of a fluid beyond the Navier-Stokes 
equations at large values of the gradients. Thirdly, they show that correlations in 
the particles positions and velocities play an important role in the physical proper- 
ties even in low density gases. The slow, non-exponential decay of the equilibrium 
time-correlation functions associated with linear transport coefficients discovered 
around 1968 with molecular dynamics simulations is one of the dramatic effect 
of correlations in dense fluids. The physical origin of those so-called long time 
tails is sometimes referred to as vortex diffusion and occurs through so-called 
mode-coupling effects. The latter occurs when a tagged particle in the fluid cre- 
ates, through collisions with and between the surrounding particles, a back flow 
of momentum through the fluid, which ultimately returns to the tagged particle. 
At liquid density, correlations lead to other dramatic effects. The concepts of 
mean-free path and of collision time cease to have a clear meaning, and the mo- 
tion of individual particles is dominated by cage diffusion: owing to pronounced 
structural correlations, each particle flnds itself trapped for some period of time 
in the cage formed by its immediate neighbors, rebounding against it and thereby 
slowing down the diffusive motion. Cage diffusion is key to understand the "uni- 
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A BBGKY hierarchy 

The BBGKY hierarchy ^ is an hierarchy of equations coupling the equal-time 
reduced distribution function /„, 

.fn{Xl, . . . , Xn ) 

= I dX„+r...dXNF{X^,...,XN;t). 

The latter can also be expression as ensemble average of microscopic fields, e.g. 

/ ^ 

MX,t) = (^<5(a;-x,(t)) 
\j=i 

/ JV AT 

h(X,X';t) = (^^5(x-x,(t))5(x'-x,(t)) 
\j=i k^j 
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and 



f3{x,x' ,x"-t) 

N N N \ 

E E E '^(^ - x,(f))5(a.' - Xfe(t))5(x" - x,(t)) ) , 
\j=l kj^i lj^j,k I 

and so forth. 

The equations of the BBGKY hierarchy can be obtained by direct integration 
of the Liouville equation The n-particle reduced distribution f„(Xi, . . . ,X„) 
represents - up to the factor A*'! / (A'' — s)! - the probabihty density of finding simul- 
taneously n particles in the specified phase-space positions Xi, . . . , X„. In practice, 
those equations are often recast in terms of the correlation functions defined from 
the cluster representation of the reduced distribution functions for n > 2, 

f2{x,x';t) = g2{x,x',t) + fi{x,t)fi{x' ,t) 
f3{x,x' ,x";t) ^ gz{x,x' ,x" ,t) + [fi{x;t)g2{x' ,x" -,1) + c.p] 
+h{x,t)h{x',t)Mx",t) 

and so on, where c.p. means cyclic permutation. In terms of the standard notations 
(llSl and (1151), the first two equations are 



dt 



Li 



/i(l) = y"d2Li2/i(l)/i(2) + Li2ff2(12) 



and 

d_ 
dt 



52(1,1') =ill'./l(l)/l(l')+ J d2[Li2/l(l)32(l',2)+ V2/l(l')52(l,2)] 

+ I d2 (Li2 + V2)/i (2)52(1,1') 

+ Ln'52(l,l') + y"d2(Li2 + Li.2) 53(1,1', 2). 



B Useful properties 

We gather here some basic properties that are often used in the remainder of the 
paper: 

— Equal-time Poisson bracket: 

[N{X,t),N{X',t)]^^ (112) 

= {S{X~X'),N{X,t)}{X). 

Therefore, 

X^/-"(X, t-X',t) = T {Six - X'), fix, t)} iX) . (113) 

— From rotational symmetry, 

j dX2Li2SiXi - Xa) = . (114) 

— Unlike their quantum counterparts, the n-point correlation functions c'"' are 
fully symmetric with respect to there arguments. 
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C First and second derivatives of f2[<p] 



In this section, we present the main steps leading to Eqs.(34| and (36 1 [?]. Since 
0[(j>] depends on through S[(j>;xQ], its functional derivatives are combinations of 
the total action's derivatives. For instance, we easily find. 



and 



50(1)50(1') 



00 



sn 



5S 



S<j}{l)S<j}{l')\4-o 

SS I 

0(1) Uo 



(115) 



(116) 



SS I 

50(1) Uo 



SS I 
50(1') l-*o 
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(i')Uo 



Wh en S is linear in the external potential, the term in 5^5/50(1)50(1') vanishes and 
Eq.(116l equals the two-point correlation of 55/50(1); this is the typical situation 
in field theory as discussed before. In the present theory, the action is non-linear in 
0, 5^5/50(1)50(1') is nonzero and generates, in addition to the correlation function, 
the term TZ{1, 1') related to the response functions x^- 

The derivatives of 5[0;2;o] at a given potential 0o can be obtained by a sys- 
tematic perturbation expansion in powers of the variations 50 around 0o. To this 
end, we introduce the following notations. With = 0o -|- 50, we write the total 
Hamiltonian as 

^0 = 'H(/)o ~t~ Sl-t , 

with SH = JdXjV{X,r)54){X,r). We denote by x{r) = (g(r),p(T)) and xo{t) = 
{qo{T),po{r)) the closed-time trajectories starting at xq and governed by the Hamil- 
tonian Hrf) and Hifig, respectively. Given a dynamical variable A, we denote by 
A — Ao + SA = Ao + Ai + A2 + ■ ■ ■ its expansion in power of 50, where An is the 
term of order n in 50. For instance, SH = SHi + SH2 + ■ . ■ with 



SHi = t 
SH2 = 



'(r) J dXMo{X,T)5<j){X,T) 
t'{r) J dXNi{X,T)S<j){X,T) 



(117) 



= qi ■ dq5Hi{qo,po,T) + pi ■ dp5Hi{qo,po,r) . 
With these notations, we rewrite the Lagrangian as 

£0(17, q, t) = [po ■ go " 'H4,{xq)\ + [6p ■ qa -po- 5q + H4,(xo) - ^0(2;)] 
+ [po ■ Sq + po- Sq] +Sp- Sq 

so that by time integration and integrations by parts of the last line, we obtain 
5[0] = 5o + ^ dr {5p ■ qo - Po ■ Sq + H4,{xo) - 'H4,{x)) 
+ - ^ drSp ■ Sq — Sp ■ Sq 



+ Po(Tf) ■ Sq{Tf) + -Sp{Tf) ■ Sq{Tf) ~ 5^(0) 



(118) 
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The boundary term is 



St[(po + S(t),xo] = Po(t/) ■ Sq{Tf) + -5p{Tf) ■ 5q{Tf) 



(119) 



so that, with (j) — 4>o + S(j>, the last hne of Eq.( 118 1 vanishes 

5S = <j> dr {6p ■ qo — Po ■ Sq — K4,{xo) — ^^(a;)) + ^ dr 5p ■ Sq — Sp ■ Sq . 

(120) 

Using the Hamilton equations for the unperturbed trajectories, the Taylor expan- 
sion of the Hamiltonian around the unperturbed trajectory gives 



'H4.0 (x) - (a;o) = -po ■ 5q + 5p ■ qo + uf^^ {5x) + , 



(121) 



and where the second order term is 



Sq 



Sp" 



00 



2 dgq'H4,g{xo) + Sq- Sp- dqpH^a(xo) + — ■ dpp'H4,„{xo) .(122) 



Using Eqs.( 121 1 and ( 122 1 in ( 120 1, the total variation of the action writes 

+ . . . 



SS 



dr 



- (5p -Sq-Sp- Sq) - H^^l {5x) - 5H{x) 



At this stage, we have removed all the terms that cancel irrespective of S(f). We now 
proceed with the expansion in orders of 5<j). To third order in 5(j), the integral in 
the last expression is 



dr 



\ (Pi -qi-pi- qi) - y-f^ixi) - SHi - S'H2 



O 



The first order variation (gi(r),pi(r)) of the trajectory satisfies 
■ gi = gi ■ dqpH^g{xo) +Pi ■ dppH^g{xo) + dpSHi{xo) 

, Pi = -gi ■ dqqH^g (Xo) - Pl ■ dpqH^g (xo) - dqSHl {xo) 



and therefore 



^ / dr [pi ■qi-qi-pi]= (j) drUflixi) + ]^<i> SH2 



(123) 



(124) 



Combining Eqs.(124l and (1231, to second-order in S<j), the total variation of the 
action is 



SS 



- j dTt'ir) j dXMo(X, t)S<}){X,t) 
-)^jdTt'{T) j dXMi{X,T)54>{X,T)^0{5(t?). 



(125) 



In order to express the second term in EqJT25|, we first note that, using the results 
of linear response theory recalled sectiontJl tlie density variation Mi can be written 



M(i) = * dV R{l,V)S(l){l') 
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in terms of the response function over the closed-time contour, 

i?(i,i') = e(^-r'){AAo(i),M5(i')}- 

Therefore 

j dlM^{l)5cj>{l) = <j> dldl'Ti{l,l')5(l>{l)S(t>{l') 



= ^jdljdl' [nil, V) + 71(1', i)] 5(l>(l')5(t,(l) 

= jdlj di'x(i, l')S<j,{l)5^{V) , (126) 



x(i,i') = ^(7;{Ar„(i),AAo(i')}. 



with 



The power expansion of the total action ( |125[ ) becomes 

s[(t>o + s<j}] = s[(t>o] - LiAAo(i)5<^(i) - I- Ildidi'xii, i')5(i>{i)s<j,{i) + 0(54^') 



2 

Therefore the first two functional derivatives of the action are 

SS 



(1) 



= -AAo(l) 



and 



5,^(1)5^(1') 



which, combined with Eqs.(115l and (116 1, yield the desired results 



D On the differentiability of the action functional. 

Although it is not at all the purpose of this paper to be concerned with math- 
ematical details, we feel that a word of caution regarding the differentiability of 
f^[(t>] and the derivatives i?'"-* is in order. Strictly speaking, the effective action 
functional f^[<j)] is not (Frechet) differentiable everywhere. It is however infin itely 



differentiable at any physical potentials <l>p. In that case, the term 56 in (29 1 can- 
cels the boundary terms resulting from the variations of i(r/) in the first term 
JJ^^ dr C^iqij), v{r), r), which allows for the definition of the functional derivatives. 
For a non-physical potential, the cancellation does not necessarily occur and the 
definition of the functional derivatives is not always possible. As a consequence, 
the quantities i?'"-* can not be regarded as functional of cj> and expressions like 
are not correct, and therefore using it could give wrong results. 
This is unfortunate since such recurrence relation is very useful to derive by simple 
differentiation the hierarchy of properties satisfied by all the successive derivatives 
from a single starting relation. 

To regain the benefits of such recurrence properties, we define the following 
formal device. Given the ensemble average {A{X, t)) of a dynamical variable A, we 
define the generalized ensemble average as 

(M(i)))0 = ^(e*^'^'^(i)) , (127) 
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for any potential ^ At 4> = 00, the generalized average reduces to the physical 
ensemble average (IM, 

{{A{X,r)))^^^, = {A{X,t{r))). 

The generalized average is a functional of <j), which, like S, is differentiable at any 
physical potential 4>v Using results of the previous sections, we find its functional 
derivative at (j>o, 



{5A{1)5M{1')) + 



<t>a 



5A{1) 



<t>0 



(128) 



in terms of the fluctuations 5 A — A— (A) . Since, using the results of linear response 
theory. 



(129) 



we find 



{%{A{l),J^(l')}) - -{SA{l)Sm')) 



(130) 



For instance, with A = A/", we shall write 

/,(i) = ((AA(i))),. 

The phase-space distribution is 



and from Eqs.( 130 1 and (|36 1. 

Ml 



(re{AA(i),AA(i')})--(<5AA(i)5AA(i')) 



= ^7<''(i,i'). 

Therefore we can obtain the successive derivatives of by writing the n-th 

derivative as 



(n) 



(131) 



in terms of the dynamical variable oj'"', and obtain the (n + l)-th derivative by 
functional derivative of its generalized average at (j) = (j)o, 



("+i)(i. 



,71+1) 



(132) 



<t>0 



This scheme can be used to systematically generate explicit expressions of the 
derivatives of the effective action [53] as co mbin ation s of p roducts and Poisson 
brackets of the fundamental field. The steps ( 131 1 and ( 132 I can be schematically 
written as 



("+!)/- 



,n + l) 



S(t>{n+l) 



-«"(i, 



with the understanding of the caveats discussed before. In this work, successive 
functional derivatives are obtained according to the scheme just defined. 
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E On the occurence of the linear combination of correlation and 
response functions 

As shown in appendix [d] for any dynamical variable A{X,t) we have 
5{A{1)} 



= I {Tc {A(l), AA(l')}) - i (M(1)5AA(1')) . (133) 

</>0 

From the expression emerges a fundamental operation between dynamical variables, 
which we find interesting to highlight. For any dynamical variables A and B defined 
on the closed-time contour, we define with square brackets [•,•] the symmetric 
operation, 

[Air), B(r')] := ^ {% {A{r), B{r')}) - J {SA{r) SB{r')) 

= [B(r'),A{r)]. (134) 

In term of this operation, equation (|133l) reads 



= [A(1),A^(1')] . (135) 

00 



and 

X(l,l')= K(1),AA(1')] . 

The operation ( |134[ ) combines two fundamental quantities of classical statistical 
mechanics, namely the Poisson bracket [■, ■]pg and the product of dynamical vari- 
ables, which are intimately related to the response and correlation functions. It 
obviously respects the dimensions thanks to the action parameter s (the Poisson 
bracket [^4, B]pb has the dimension of the product AB divided by the product of 
a length times a momentum, i.e. by an action.) The Poisson bracket carries the 
dynamical information on how fluctuations are propagated in time, while the cor- 
relation function gives information about the likelihood of fluctuations (occupation 
number.) At equilibrium, both quantities are simply related by the fluctuation- 
dissipation theorem and only one is useful (one usually works with the phase-space 
density correlation functions.) The emergence of the linear combination provides a 
justification of the fact recognized by MSR and others [H that a time-dependent 
theory of classical dynamics should be expanded to involve not only correlation 
functions but also response functions in order to be amenable to a systematic 
renor maliz at ion . 

We can shed some light on the occurrence of that particular combination by 
regarding the action function as the classical limit discussed above (see discussion 



after Eq.(31l.) In a quantum system, both fluctuation statistics (correlation) and 
response effects can be characterized by products of the same quantum field opera- 
tor; the intensity of fluctuations involve the anti-commutator [•,■]+ and the response 
is related to the commutator [•,■]-. In practice, both properties are contained in 
the time-ordered Green's functions, say G(r, r') = {T^P'(r)!Z'(r')) , where !f'(r) is 
the field operator (for simplicity, we do not write the other degrees of freedom). 
We can re-express the Green's function in terms of the averaged commutator and 
anti-commutator of the field, 

G(r,r') = {i{T)^{T'))e{T - r') + {i{T')^{T))e{r' - r) 

= ^(^[ir(r)Mr')]_) sign(r - r') + ^ ^ [#(r), *'(r')] . 
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1' 2' 




V 2' 



(5"+" n 



/.p(l)...%(n)#A(l')---^<^A(n') 



= f^i...i2...2(i---";i'---«') 



Sfp{l)...Sfp{n)5Ml')...SU{n') 



= ri...i2...2(i---";i'---«') 



Fig. 3 Representation of derivatives of the effective effective action f2 and effective 
potential F. An inward (outward) arrow indicates derivative with respect to the 
physical (non-physical) component. 



From the correspondence relations, in the classical limit, the anti-commutator [■,■] + 
becomes a product, the commutator becomes the Poisson bracket ift{ip'(r), !i'(r')} 
and therefore 

G(r,r') ~> ^ (r [l'ir),1'(r')]^^) + (^(r)tf'(r')) 



where we recognize the bracket |l34| introduced earlier. 



F Closed vs unclosed expressions for ^(1, 1') 



The previous closed expression for |x'(i, 1'), /(i')|, Eq.(60l, should be compared 
to its "unclosed" v ersi on, which can straightforwardly obtained by using the equa- 



tions of evolution ( 47 1 and ( 49 1 for x s-nd J?'-"'' on both sides of the equation ( 57 1 
for S. We obtain, 

{r(i,i'),/(i')}(i') 

= ^73(123)53(231') (136a) 

+ (^|)'73(i23)73(i'2'3')[x'''(232'3') - ^x(22')x(33') - x'''(234)x-'(45)x<''(52'3')] 

(136b) 



+ Zi(r')i73(i23)x(2,2)x(3,3)^ [c(i', 2')x-'(2', 2) 



S_ 



(136c) 



thus establishing the correspondance between Eqs.( |60b| ) and ( 136b[ ) 

In his "fully renormalized" kinetic theory of equilibrium liquids, Mazenko [12] 
uses as a starting point an e xpression for the memory function equivalent to the 
equilibrium limit of Eq.( 136 1. 
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1' 2' 



5f,{n + 1) 



1' 2' 



1' 2' 



Fig. 4 Graphical representation of chain rule. 



G Derivation of Eq.{92) 



Equation (91 1, namely 



{r^(i, 1'), /(I)} (1') = f (r> (t, t') + r< (t, t')) + ix^,x^;t)5{t' - to) 

+<5(t-to)C(l,2)r^(2,l') 
implies that F'^ is of the form 

1') = r'-^«(i, 1') + Sit' - io)r*^(Xi, x(; t) + sit - t„)r^'\x^,x[;t') (137) 



where J^''^^ and are given by Eqs.(93l and F'^' is a regular function of time 
variables statisfying, 



{r^^HXi,X[;t'),f{X^,to)} = C(Xi,to;2)r^(2,l'). 



(138) 



_r'^' can b e rela ted to and the initial correlation function Co as follows. Intro- 
ducing Eq.||137|l into (IstI yields 



C(l,l')= dt3 // dX2dX3x''(l,2)r'-=^(23)x''(3,l') 

J to J to J J 

+ f^dt2jj dX2dX3x-"(l,2)r*(X2,X3;t2)x'^(^3,to;l') 

+ JdXiJ d3x''(l;X2,to)r(^>(X2,X3;t3)x'*(3;l') 

Setting t' = to in th e pre vious expression, the resulting expression is used to cal- 
culate the rhs of Eq.( 138 1, which gives 



{F^'\X^,X[;t'),f{X^,to)} 

= {F'{X[,Xv,t'),fil, to)} + S{t' - to) J d2x''(Xi, to; 2)f'^^^ {X[, X2, t^) . 

Using x^('f 1, io; ^2, to) ~ \^f(Xi),S{X\ — X2)}, the previous relation implies the 
following expression for _r'^' 



F'^^\XuX[-t')^F'{X[,Xr,t')+5{t' ~ to)r(o' [X^ ,X[), 



(139) 
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3 = 5(1-3) 
3 = 5(1-3) 



3=0 



Fig. 5 Graphical representations of Eqs.( 142a I, (142b I and (142c 



with r^''^ given by Eq. ( 94 1 . Combining Eqs. ( 139 1 and ( 137 1 yield the desired result 
Eq.(p2l. 



H Diagrammatic representation 



H.l Derivation of Eq.(104) 



From the chain rule of differentiation, we have 

1^ =5(1-3)= / d2ri2(12)f2i2(23) + rii(12)r222(23) 
-0 = y d2r22(12)«i2(23) + r2i(12)l222(23) 



5/p(l) 



and 



=0 = / d2Ai(12)r22i(23) 
^yfll = 5(1 - 3) = I d2r2i(12)r2i2(21) 
where we used =0. We deduce Ai =0, and therefore 



5(1-3) = J d2ri2(12)J7i2(23) 
5(1-3) = J d2r2i(12)J7i2(21) 

J d2r22(12)J7i2(23) + r2i(12)r222(23) 



= 



Combining equations (142 1 with (102 I leads to Eqs. (104 1 



(140) 
(141) 

(142a) 
(142b) 
(142c) 



H.2 Higher-order derivatives with diagrams 

We define here an economical graphical representation for deriving relations be- 
tween the correlation and response functions with the vertex functions. We define 
the graphical notations as shown in Fig.([3|. 

Using the chain rule, the functional derivative S/Sfp and S/SJa of 
add an additional two-point vertex_to the diagram as illustrated in Figs|4] and js] 

vatives by 



Using Ai = and the rules in Fig.(|4|, we can evaluate the successive derivatives ay 
letting 5/5fp^A act of the diagrams of Fig.(|5|. We illustrate the calculations leading 
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Fig. 6 S/Sfp of the first diagram in Fig. ([5 





Fig. 7 5/5fp of last diagram in Fig.([& 



to tlie tiiird-order derivatives Eq.( |105a| used in the body of the paper. S/Sfp of the 
first diagram in Fig.([5| yields the diagram in Fig.S, and S/Sfp of last diagram in 
Fig.|5| yields the diagram in Fig.([7]l Combini ng the diagrams in Figs. ||6| and ^ 
yields the diagrammatic representation of Eq.(105al shown in Fig.®. 



I Details on the derivation of the real-time formalism 

By considering the equation of motion ( |43[ ) for Af{r) in an external potential (f> on 
both sides on the closed-time consider, it is straigthforward to extract the evolution 
equation for the classical and non-classical components Afp,A as 



LI 



Np{l)-\LtNA{l) = l73(l,2,3)(7V,(2)iV,(3) + ^iV4(2)iV^(3)) 



LI 



iV4(l) - iVj,(l) = 73(l,2,3)iVj,(2)iV4(3) 



with the single particle Liouville operators, 



L?. = {/io(l) + <^p(l),.} 
Lt» = {0^(1), •} 
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\ (2) (2) CO) 

/ @ @ @ 



Fig. 8 Figs.(|6|) and (|7|) can be combined to lead to this diagram 



By averaging these equations over the initial conditions, we obtain the evolution 
equation of for the averaged densities fp^A = {J^p,a), 



dti 



LI 



(143) 



^73(1,2,3)(/p(2)/p(3) + ^/a(2)/4(3)) + |73(123)r222(23) + 5(t-to)/p(l) 



and 



dti 



LI 



fA{l) - Lffp{l) = 73(l,2,3)/p(2)/zi(3)is73(123)fi2i(23) . (144) 



As in Sec. |3.6[ equations of evolutions for the second-order derivatives Qij, i.e. for C 
and x^'^ y can be obtained by differentiation of Eqs.( "l43|l44 1. The terms involving 
the third-order derivatives Oijk in the resulting equatioiis (not shown here) can be 
expressed in terms of memory function kernels X'^''*''-^ such as 



is ( 2^212(231') 

2 73li23j I ^^^^^231') 2^2222(231') 





-02 



(1, 



^12, , 



which lead s to eq uations (106 1 and (1071. Following the same steps as in Sec. 3.7 
using Eqs.( |85]|86t , the evolution equations for x^''^ 8-nd C yield the following (non- 
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trivial) relations between E^''^''^ and p^-^''^ , 
{r2n(123),/(l)} = {ri2(32),5(l ^ 2)} - {r2i(12),5(l - 3)} - 7123 
{ri2i(123), /(I)} = {ri2(32), 5(1 - 2)} - {ri2(12), 5(1 - 3)} - 7123 



5/p(3) 



{r22i(123),/(l)} = -{r22(12),5(l-3)} 



1 5X''^(12) 
is Sfp{3) 



-5(ti-to)^^[r222(14)A2(42)] 
{r222(123), /(I)} = ^{r2i(32), 5(1 - 2)} - 1{ A2(12), 5(1 - 3)} - ^7123 



1 5r^(i2) ^^^^ _ t^)_A_[a,,{u)ri2m] 



is 5/zi(3) 



5/4(3) 



{ri22(123), /(I)} = {r22(32), 5(1 - 2)} - {r22(12), 5(1 - 3)} - 7123 



5^-^(12) 



{r2i2(123),/(l)} = {r22(32),5(l-2)}-7i23 



5r^(12) 
5/4(3) 



Those relations are then combine d wit h EqJ K)7l to close the hierarchy through the 
functional differential equations ( 108 and irairor the memory function kernels. 



